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ABSTRACT 


This thesis is addressed to the development of fast and 
efficient methods for voltage contingency selection. One of the 
problema faced in the real time execution of security analyais^ 
especially the voltage security analysis, is the non-availability 
of fast and accurate methods for predicting the post outage 
conditions. Hence, the attempts have been made in this thesis to 
suggest a new set of voltage and reactive power distribution 
factors which can be used for rapid coaiputation of' the post outage 
system states with sufficient accuracy. 


Another difficulty is due to non-availability of proper- 
performance indices which will reflect the true relative 
severities of cont ingenci es . Various contingency ranking methods 
available in literature, in general, suffer from the masking and 
mlsranking effects. To overcome these problems, several new 
higher order performance Indices have been explored. A simple 
approach, based on least square error minimizat iot^ to compute the 
optinial values of weights associated with voltage and reactive 
power performance indices has been suggested. 


The 

model and 
sel ect ion 


potential of the proposed distribution factor based 
the new performance indices for voltage contingency 
have been tested on IEEE-14 bus and a 19 bus Indian 


system. 



-o. . . j J 4(i.Vi) 

£ (‘^<^2.- “ R.)A^0 ^ j" 




DEDICATED TO 


MY PARENTS AND AUNT 



ACKNOWLEDGEMENTS 


With profound sense of gratitude I express my sincere thanks 
gurus and guides Dr . S. C, Srivastava and Dr.P.K.Kalra for suggestin 
topic and guiding' me throughout course of the thesis work. My 
have be^en the constant source of inspiration and encourageme 
times when I was passing through difficult momeiits.My special 
are due to Dr . S. C. Srivastava for sparing his valuable time to com 
the thesis work in--time. 


1 am grateful to Dr . L . P . Singh , Dx- . R . Ai'ora , Dr . S . S . 

Dr . R . K. . Varma , and Dr . P. R. K. Rao for teacdiing me various adv 
courses dui.'ing my M.Tech, I would like to express my hear 
thank.s to Mrs ( Dr.) Damayaiiti Singh for teaching me the adv 
features of Ehiglish language. 

I am fortunate to be in the splendid company of Mahendra, 
Frasad, and Chandra Shekax' in. particular and all H-top Junt 
general# during my stay here. My alter ego Ved needs a sp 
mention for he being my Hindi guru. 

My special thanks are due to K . N . Srivastas'^a for the val 
discussions we had during our course work and thesis work.: 
grateful to S.N. Singh for pt'oviding me UPSEB data and it’s det 
I am thcinkful to my Doordarshan friends, Defence frieiK-l.s and a: 
my friends who involved directly or indirectly in thesis px'epar 
and gave me good company in the campus . 

I am indebted to Capt. Rajeev Sawhney and his family 
•Dr . N . N . Kishore and his family who made my s tay at IIT K memox'able 

I am also thankful to Mr. Bajpai for his fast and efficient 
writing. 


’ - R;3io 



CONTENTS 


Page Mo« 

ABSTRACT 

CHAPTER ” 1 INTRODUCTION 1 

1.1 General Introduction 1 

1.2 State of the Art and Motivation 6 

1.3 Theaia Organiaation 13 

CHAPTER - 2 VOLTAGE AND REACTIVE POUER DISTRIBUTION 14 

FACTORS 

2.1 Introduction 14 

2.2 Reactive Power Flow Model of Tranamisaion 16 

Uni t a 

2 ; 3 Proposed Distribut ion Factors 19 

2.4 A Sensitivity Baaed Approach to Compute 22 

Dia trlbut ion Factors 

2.4.1 Line Outage Voltage Distribution 23 

Factors 

2.4.2 Generator Outage Voltage Diatri- 26 

bution Factors 

2.4.3 Line Outage Reactive Power 28 

Distribution Factors 

2.4.4 Generator Outage Reactive Power 29 

Distribution Factors 

2.5 System Studies and Results 30 

2-6 Conclusions 32 

CHAPTER - 3 PERFORMANCE INDICES FOR VOLTAGE 42 

CONTINGENCY SELECTION 

3.1 Introduction 42 

3.2 Existing Performance Indices 44 



Page No. 


3.3 Pcopoaed Performance Indices for Voltage SO 

Contingency Ranking 

3.4 Proposed Scheme for Ueights Adjustment 52 

3.5 System Studies and Results 55 

3.6 Conclusion 5S 

CHAPTER - 4 CONCLUSIONS 69 

REFERENCES 73 

APPENDICES 77 

APPENDIX - A , Data of IEEE-14 Bus System 77 

APPENDIX - B Data of a Practical 19-Bus (UPSEB 400 kV 
Network) Indian System. 


80 



CHAPTER - 1 


INTRODUCTION 


1.1 GENERAL INTRODUCTION 

Security is a term used, to reflect a power system's ability 
to meet its loads without unduly stressing its apparatus or 
allowing network variables to stray from prescribed ranges. Of 
particular concern are the sudden changes or disturbances in the 
system. The framework normally used to examine the security of 
the power system network is baaed on its probable transition of 
operating states following a disturbance or contingency. 

The power system is thought of as being run under two seta 
of constraints -. load and operating. The load constraints impose 
the requirements that all loads be met, the operating constraints 
impose lower and upper limits on network variables. The system 
states have been broadly classified as Normal, Emergency and 
Restorative states. The system is in the 'normal state' when the 
lead and operating constraints are met. The system is in 
'emergency state’ when there are violations of operating 
constraints, though satisfying load constraints. In 'restorative 
state' the system operating constraints are met but load 
constraints are not satisfied. Further detailed classification of 
states are reported in literature such as a five state model of 
Fink and Carlson [5]. 
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Security la defined with, respect to a set of random events 

called the 'aet-of-next-contingeneies ' . This is collection of 

disturbances that could happen. A system is said to be secure if 

it is in the normal state and none of the contingencies could 

cause its transition to the emergency state. If all the possible 

contingencies in a power system network (y.tB considered, it can 

never achieve a 'secure' state. Hence, the security is referred 

only with respect to a certain preapeclfied credible contingencies 

whose probability of occurrence la high. These disturbances, are 

either in the form of 'network outages' such as a line or a 

transformer outage or in the form of 'power outages' e.g. a 

generator outage. Following these outages, the system becomes 

deficient of transmission and/or generation capabilities, causing 

the lines to be overloaded beyond its thermal or stability limit 

and bus voltages. These two violations viz. line flows and bus 

A 

voltages are mainly related to the deficiency of real power flow 
and reactive power support respectively. Hence the two are 

studied as separated subproblems known as 'Line flow or MU 
Security' and 'Voltage security' studies. This thesis mainly 

concerns with the 'Voltage security' problem which has attracted 
the attention of researchers in the recent past. 

Secure operation of power system requires the assessment of 
security with respect to the preselected set of contingencies and 
the planning of corrective control actions if it is found to be in 
'insecure' state. The two major functions of power system 
security are i) Security assessment and li) Security control. 
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Security analysis and control have been implemented through 
a number of software packages in modern energy control centers. 
The major components of on-line security analysis are shown in 
Fig. 1.1. The monitoring component starts with the real-time 
measurements of physical quantities such as line power flows, line 
current flows, power injections, and bus voltage magnitudes; as 
well as, the status of breakers and switches. The measurement 
data are telemetered from various locations to the control center 
computer. Bad measurement data are rejected by filtering the 
transmitted data through a simple check of their reasonability and 
consistency. The remaining data are first systematically 
processed to determine the system configuration (generator and 
transmission network connections) or network topology. Then the 
available data are further processed to obtain an estimate of the 
system state variables (bus voltage magnitudes and phase angle for 
normal steady-state). State estimation is a mathematical 
procedure for computing the ’’best” estimate of the state variables 
of the power system baaed on the available data, which are in 
general corrupted with errors. 

To assess whether a normal operating state is secure or not, 
a set of contingencies are analyzed. Therefore, to assess the 
system response to contingencies, a contingency evaluation is 
carried out using the on-line load flows. The on-line load flow 
uses the actual load flow model of one’s own system (from the 
state estimation solution) together with a system representation 
of the unmonitored network and neighbouring systems, i.e., an 




Secure Insecure 

slate stole 


Fig-.l-l Security monitoring and analysis functions 










external network model. Because the contingencies are future 
events, a bus load forecast is needed. 

In a practical power system, the number of contingencies to 
be considered are so large that it is not possible to analyze all 
of them in real time requirement using on-line load flows based on 
full AC load flow models. Hence all the contingencies are first 
analyzed using some approximate load flow models or using 
precomputed distribution factors and ranked according to their 
severity. In order to assess the relative severity of the 
contingencies, scalar performance index is calculated for each of 
them which is function of the variables being monitored such as 
line flows in case of 'MU security’ and bus voltage or reactive 
power injection in the case of 'Voltage Security'. The process of 
analyzing the contingencies using approximate models and ranking 
them according to their relative severities is known as 
’’Contingency Selection”. Once the contingencies are ranked 
according to their severit^*iij , the contingency analysis is 
performed using full AC load flow model starting from the most 
severe case, and stopping at a contingency which does not cause 
violation of operating constraints. Uith this two step procedure, 
normally a f ew contingencies are required to be analyzed by full 
AC load flows, thus reducing the time for security analysis. 
However, on-line security analysis has remained still a challenge 
to the power system researchers. Some of the desired features of 
a contingency selection algorithm are : 
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i) It should employ an extremely fast and sufficiently accurate 
method to compute the post outage state of the network for 
each of the contingencies. 

ii) It should employ a ranking or screening method which 
accurately reflects the severity of the contingency cases. 

1 . 2 STATE OF THE ART AND MOTIVATION 

Power system security, being an important aspect of modern 
day power systems to achieve a high degree of reliable electric 
supply, has attracted researchers since past two decades and has 
created renewed interest in the recent past. 

It la not possible to review the large number of literature 

available in the area of power system security. Hence a limited 

literature survey of the works related to the topic of the present 

thesis has only been carried out. Moreover, a recent paper by 
Balu et al . [30] provides the state of art in the area of power 

system security analysis and presents an overview. 

In a practical large power system there can be thousands of 
credible contingencies. Keeping in view the enormous 

computational effort involved, the total number of contingencies 
must be reduced for detailed investigations. This is indeed the 
purpose of contingency selection. Contingency selection itself 
requires fast computational approach because of the need to study 
thousands of cases in a short span of time. The severity of an 

outage is indicated by the numerical value of the performance 

index PI, which is a function of the line overloads or bus voltage 
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deviations. Hence, linear, non-iterative approximate techniques 
have to be employed for such real time applications. Inverse 
matrix modification Lemma (IMML) [4] and compensation techniques 
have been presented in the literature for this purpose. For line 
or MU contingency ranking Generation shift distribution factors 
(GSDFs) and line outage distribution factors [1] are very popular 
because of their speed of solution. 

G.C. Ejebe et al . have presented a fast technique for the 
automatic ranking and selection of contingency cases for a power 
system contingency analysis study [6], A contingency list is 
built containing line and generator outages which are ranked 
according to their expected severity as reflected* in voltage level 
degradation and circuit overloads. Besides this paper, many 
efficient and reliable algorithms were developed [8,9,13,21] . But 
most of these techniques can only be applied to MU limit security 
problems. 

On tb«f other hand, voltage problems also forms a very 
important aspect of security assessment. Many papers have been 
published in literature on voltage security assessment [7,10,12, 
14,15,17,19] . 

If a contingency causes voltage problems, then, in may 
instances, these problems are in the neighbourhood of the dropped 
element- Utill^iing this property, Zaborsky et al . developed a 
method [7] which uses a Gauss-seidel solution algorithm, starting 
the solution on buses near the outage, and expanding outward. 
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During the solution, voltage on buses near the outage are updated 
while the voltage on buses distant from the outage are assumed to 
be constant. As the solution continues, buses farther and farther 
from the outage have their voltages, until a point is reached 
where the changes are insignif icant . Voltages of buses still 
farther from the contingency can be assumed to have remained 
constant. The drawback with this method is that the algorithm may 
overlook certain changes in bus voltages remote from the 
cont ingency . 

Khu et al - [12] have reported a fast non-iterative 
linear izat ion method devised to evaluate the effects of single or 
mul t ipl e-cir cui t cont ingenci es upon the system laod-bus voltage. 
The standard DC load flow technique has been used to estimate the 
change in load-bus voltages due to the outage of one or more 
circuits. The authors also presented a method to calculate the 
reactive power change required to maintain scheduled voltages at 
generator buses. The method assumes as many as eight assumptions 
and some of them are not justified in a true power system network. 

Lauby et al . [15] have compared the three methods of 
contingency selection - a performance index method, a local 
solution method, and a single iteration of the Stott-Alsac 
decoupled load flow. The authors observed that for large number 
of contingencies, the PI method which is a function of real flow 
in lines would be the fastest of the ranking methods. The method 
performs well in detecting widespread and severe voltage problems. 
However, the PI method is not well suited for the detection of 
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local problems. Overall the local solution method has been found 
to be well for selecting circuit outages. In fact it performed 
better than the single iteration of the decoupled load flow. 

K. Nara et al . [17] have presented a ' new concept in 
formulating a performance index for contingency selection 
concerning voltage security analysis. The proposed performance 
index has the ability to select the severity of voltage limit 
violations. The infinite norm is used as a filtering algorithm to 
rank the contingencies in a paper published by Uasley et al . [14]. 
A method for evaluating the effectiveness of the presently 
automatic contingency algorithms in capturing contingencies which 
give out-of-limit conditions has been presented by Halpin et al. 
[16]. 

Stott et al . [18] presented various enhancements to the 
automatic contingency selection approach have shown how general 
ranking formulae can be extended to the multiple branch 
contingencies . 

First time in the literature, Marija Illic-spong et al . have 
tried to develop a new formulation of distribution factors which 
is suitable for the analysis of the reactive power problem [20]. 
They used the S-E graph and its decoupled version Q-V graph to 
define the new distribution factors. Followed by their work 
Taylor et al . [28] have presented an algorithm that can be used in 
analysing reactive power flow contingencies. This approach uses 


the widely used MU distribution factors in conjunction with 


10 


another set of VAR distribution factors to solve interatively for 
the post-contingency bus voltage magnitude changes. But these 
algorithms [20, 2S] utilize iterative scheme and do not provide a 
set of distribution factors which can be used for direct 
calculation of post outage states. Recently, Lee and Chen [29] 
have presented a method to calculate set of distribution factors 
of reactive power flow for transmission line and transformer 
outage studies. The derivation of these distribution factors are 
based on fast decoupled load flow equations. 

A fast security analysis technique for voltage security 
assessment has been presented in ref. [22] The method identifies 
the location of buses with potential voltage problems and thereby 
defines a voltage sensitive subnetwork for contingency screening. 
This allows the evaluation of large number of contingencies. V. 
Brandwajn et al . [23] have developed an efficient . contingency 

screening method for detecting both branch MU flow violations and 
bus voltage limit violations. The efficiency of the method has 
been derived from bounding criterion which reduces the 1) number 
of branch-flow computations and limit checking, and 2) number of 
buses for which the has to be calculated. The method 

is based upon the incremental angle criterion and the fast 
decoupled power flow. 

A direct ranking algorithm for contingency selection taking 
into consideration voltage security problems has been presented by 
Chen and Bose [24]. In direct ranking, the second order PI for a | 
contingency case is formulated in such a way that detailed 



knowledge of post-contingency voltages at each individual bus is 
not necessary. 

Schafer et al . [25] have presented a systesaatic analysis of 
masking effect in contingency selection and described a new 
concept with the aim to compensate for it. In this paper, the 
authors have introduced a vector norm based PI formulation, which 
allows a quantitative description of the masking effect. In ref. 
[25] a new algorithm has been presented for security analysis 
considering voltage security problems. The approch is baaed on a 
newly developed adaptive pre-filter. The authors used reduced 
models for both contingency selection and analysis. 

Schafer et al . [27] have presented an expert system based 
novel approach for contingency selection [27]. To avoid the 
masking effect which is common in algorithmic contingency 
selection methods the expert system CONTEX incorporates operator 
knowledge into the contingency selection process. 

From the limited literature survey presented above it is 
found that most of the papers on voltage contingency selection 
utilize iterative schemes or decoupled load flow models. The 
distribution factors have been popularly used for real power or 
line security problem, which is computationally much faster than 
other models. However, no such factors are available for the 
voltage security problem except the one recently suggested by Lee 
et al . [29]. The factors suggested are baaed on fast decoupled 
load flow approximattlons and thus provide quite inaccurate 
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results. Moreover the factors have been suggested to compute only 
voltages in case of line contingencies. 

Further, almost all works on the ranking methods employ the 
second order performance indices which suffer from the 
'miaranking' and ’masking effects'. Masking effect is referred to 
declaring a contingency case, having few operating constraints 
violating their limits by large amount, less severe to the cases 
having more number of operating constrainsts violating by much 
leas amount or having insignificant violation. Misrankirig of 
contingencies ■ are mainly due to the Inaccuracies in the model used 
for predicting the post outage state (voltages) of tiie system. 

Hence, the motivations behind the works carried out in this 
thesis are : 

i) To develop a new set of distribution factors to directly 
compute post outage voltages and reactive power outputs of 
sources following a line or generator outage. The factors 
are derived from the network sensitivity properties of 
Newton-Raphson method Jacobian found at the end of a base 
load flow results. 

il) To explore new higher order voltage and reactive power 
performance indices which would eliminate the masking effect 
and a new method for optimal adjustment of weights used in 
the performance indices to eliminate the misranking problem. 



13 


The ef iveneaa of the proposed distribution factors and 
performance indices for voltage contingency selection have been 
demonstrated on two sample systems. 

1«3 THESIS ORGANISATIOM 

The thesis consists of total f our chapters . 

The present chapter introduces the power system security and 
voltage contingency problem, presents a brief state of the art 
aiainly on voltage contingency selection problem and sets the 
motivation behind the present work. 

Chapter 2 presents the development of the new set of 
distribution factors which can be used for computing the post 
outage voltages and reactive powers following the outage of a 
1 ine/ transf ormer or a generator. 

In Chapter 3 new high order voltage and reactive power 
performance indices have been defined keeping in view to reduce 
the masking efect. A method based on least square error 
minlmiisat ion has been developed to compute optimal weights of the 
performance indices which reduces the misranking effect. 

Chapter 4 covers the main conclusions of the thesis and 
suggestions for future scope of research In the area of voltage 
security contingency selection and analysis. 



CHAPTER - 2 


VOLTAGE AND REACTIVE POWER DISTRIBUTION FACTORS 

2.1 I NTRODUCTI ON 

Security analyaia for real time monitoring and control of a 
power system is still a challenging task to power system 
engineers .. Outage of a transmission unit (line or transformer) or 
a generator may lead to overloads in other healthy branches and 
generators and/or cause sudden change in the system bus voltages. 
Operating personnels must know which of the outages will cause 
flows or voltages to fall outside limits in order to plan for 
preventive actions. Therefore, fast and accurate methods are 
necessary to predict the post outage effects in the system. For 

analysing real power security several fast methods have been 
suggested in the literature. Amongst them the one based on 

decoupled load flows [4] or linearized load flow such as [3] and 
distribution factor derived from D.C. load flow equations [1] are 
popularly used. The distribution factor methods and D.C. load 
flows have been widely used in real-time operation and system 
planning studies to identify branch overloads mainly because of 
the post outage state of the system can be calculated extremely 
fast and the approach la simple and direct. In many cases only 
real power flows may not be adequate in assessing power system 
security and contingency analysis, in addition, voltage magnitude 
also becomes a critical factor. Of late the problems associated 
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with r<&a.ctive power flows and bus voltages have acquired greater 
importance. The transmission capacity may sometimes be limited by 
reactive power conaiderat ions . In a few instances voltage 
collapse in a transmission network has been attributed to abnormal 
reactive power flow patterns. 

For voltage and reactive power contingency analysis, various 

¥ 

methods based on A.C. load flow have been repOted in literature. 

A 

Full A.C. load flow methods such as Newt on-Raphson load flow 
(NRLF) or Fast Decoupled Load Flow (FDLF) can not be used for real 
tittle applications because of the large computat ional time. Some 
of the approxifliat e and fast methods include the linearized load 
flows, using only one iteration of NRLF or FDLF. These methods, 
in general provide quite inaccurate results. Some of the attempts 
to derive the distribution factors for the reactive power or 
voltage contingency analysis include the efforts of Illic Spong 
and Phadke [20] and Taylor et al . [2S]. Their models are based on 

decoupled Q~V equations which computes post outage voltages and 
reactive power using an iterative scheme. They, however, did not 
suggest a set of factors which can be used for direct calculation 
of post outage state. Moreover the inaccuracies in predicting 
the voltage is as high as 10-30%. Thus, there seems to be no 
proper attempt to derive set of distribution factors for reactive 
power and voltage contingency except a recent attempt by Lee and 

Chen [29]. They suggested only voltage distribution factors for 

line outages derived from FDLF equations. Since FDLF model itself 
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ivolves several approximations, the factors suggested are bound to 
provide quite inaccurate results. 

Hence, in this chapter a new set of distribution factors 
which can be used for direct computation of bus voltages as well 
as reactive power outputs of sources following a line/transformer 
or generator outage have been suggested. The voltage and reactive 
power distribution factors have been derived using base load flow 
results and an efficient method exploiting the sensitivity 
properties of NRLF Jacobian available at the end of base load 
flow. The accuracy of post outage results using the proposed 
distribution factors have been established on IEEE-14 bus and a 
practical Indian system. 

a. a REACTIVE POWER FLOW MODEL OF TRANSMISSION UNIJS 

Two main components of transmission units are transmission 
line and transformer with on line tap changing (OLTC) provision. 
Outage of any of these units causes redistribution of reactive 
power flow in the network and hence the changes in the voltage 
profile. In order to arrive at a direct method for calculation of 
new voltage profile and reactive power output of sources, the 
effect of change in reactive power flow of tranamission unit 
considered for outage from its previous value to zero is required 
to be studied. Since the reactive power flow in these units at 
receiving and sending end are different, neither of them can be 
used to define the outage distribution factors. Hence the average 
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or t ranainl 1 1 ed powers have been referred for this purpose as 
defined below. 

Consider rr- equi val ent model of a transmission line-l 

connected between bus i and j as shown in Fig. 2.1(a). The 

complex power from bus i to j (S, . = P. , + 1 Q, ,) can be 

ij -'^ij'^ 

expressed as. 



ij 


j °ij 


= V. 
1 


[K - "iJ * ''i [■' ®cap)] 


where Y. . = (G. . j B ) ia the eeriea admittance and J B is 
iJ — cap 

the half line charg^ing admittance of the line. Equating the 
imaginary parta in the above equation, the expression for reactive 
power flow can be written as, 


“ij 


ffi. . + B Iv? - V.V. G. . sinfd.-di.l + V.V. B.. cosfdf.-d.l 
L iJ capj 1 ijij lijJ iJJ-J l^ jJ 


(2.1) 

Similarly, the reactive power flow from bus-j to bus-i can be 
expressed as. 


Q . . = 


- Tb. . + B Iv^ - V,V. G. . sinfe.-©,! + V.V. B,, cosfei -d.l 
I ij capJ J i J ij I J ij 1 J ij L j ij 


or 


Q., = - (b, . + B Iv^ + V.V. G. . sinfd.-d.l + V.V. B.. cosffif.-e.l 

“ji t°ij capj j i J ij L i jj i J ij L i jj 


( 2 . 2 ) 



u 



Fig. 2.1a Transmission line n- equivalent model. 


j 



B'ig. 2.1b OLTC Transmission n-equivalent model. 



The average tranainltting reactive power flow in the above line 
I, can be defined as 
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2 




Frofli Eqna . (2.1) and (2.2) 



1 

2 


I 1 j capj j 



V.V. G. . sin 
1 J 1 J 




(2.3) 


Fig. 2.1(b) shows the n-equival ent model of a transformer 
with OLTC provisions; 'a' is the tap setting value and ^ ^ is the 
aeries admittance of the transformer. Taps are assumed to be 
provided between transformer and bus-l. The average transmitting 
reactive power in the transformer having equivalent rr-model as 
shown in Fig. 2.1(b) can also be derived similar to the 
transmission line. 


2.3 PROPOSED DISTTRIBUTION FACTORS! 

In the voltage security or reactive power security studies, 
the effect of 1 ine/transf ormer or generator outages on the system 
voltage profile and reactive power output of sources are required 
to be studied. The outages can be either of a single unit at a 
time or multiple units simultaneously. However, in the present 
study, only single unit (line or transformer or generator) outage, 
which is the moat common, have been considered. Since the 
modelling of transformers are similar to the line, their outages 
also have been referred to as line outages. The set of proposed 
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distribution factors for line and generator outages are defined as 
follows . 


Consider the outage of a line/transf ormer - I carrying the 

average or transmitted reactive power Q^. The factor to compute 

change in voltage magnitude (say AV. at a bus-i) has been termed 

1 

as Line Outage Voltage Distribution Factor (LOVDF) a, . and is 
defined as 


AV. 

a 1 ” 1 f ■ N and i ” 1 ^ ^ ^ 


C2.4) 


For the outage of above line-l for computing the change in the 
reactive power output of source-j (AQG^) is another set of 
distribution factors termed as Line Outage Reactive Power 
Distribution Factors (LOQDFs) are deined as 

AQG . 

C, . = j = 1 , . . ,N_ and 1=1, ..,N (2.5) 

■’ Ql ^ . 


where N, Nj^ and Ng are number of buses, lines and reactive power 
sources (including generators, & synchronous condensers) in the 
system. 


Now instead of line outage, consider the outage of a 

generator - g carrying a power QG during pre-outage condition. 

& 

Due to this outage assume that the change in voltage at bus-1 is 
AV^ and the reactive power of source-j is AQGj (j f>= g) . These 
changes in voltage and reactive powers can be computed with the 
help of distribution factors termed as Generator Outage Voltage 
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D i 8 1 r i bu t i oji Factor (GOVDF) and Generator Outage Reactive 

Power D i a 1 1* 1 bit t i on Factor (GOQDF) d . respectively, defined as 

d J 

f ol Iowa 


AV. 



^ 

i = 1, . . ,N 


S = 1, . - ,N^ 

(2.6) 

and 




AQG 

“sJ ' 557" ^ 

o 

& (2.7) 

6 — 1 , - . , Nq 

The above four aeta of distribution factors can be computed 
directly from base case load flow results utilizing sensitivity 
properties of Jacobian matrix as described in the subsequent 
section. These factors can be directly used to compute the post 
contingency system voltage profile and reactive power outputs of 
sources. If the post contingency voltage at bus-i is and the 

reactive power output of source-j is QG^^ following outage of 
line-l or generator~g, and their base values are V? and QG? 
respectively, then the post outage values using the distribution 
factors can be computed as, 

For o triage 

V? = V? + a, . . oj 1 = 1, . . ,N (2.8) 

i 1 It I 


= QG^ 




(2.9) 
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For Gen-g outage 




’'el 



1 = 1 , . . ,Kr 


( 2 . 10 ) 



QG 


+ d 


e.j 


OG. 


& 


j = 1, - - ,Ng 

^ £ 


( 2 . 11 ) 


2.4 A SENSITIVITY BASED APPROACH TO COMPUTE DISTRIBUTION FACTORS 


For computlne voltage and reactive power distribution 
factors as defined in Eqns . (2.4) to (2.7), the base case line 
flows, output of sources and the effect of each outage on bus 
voltage and reactive power outputs of sources are required to be 
known. A new approach to compute the post outage changes in bus 
voltage magnitudes and reactive power generations, utilizing the 
sensitivity properties of the Newton-Raphson base case load flow 
Jacobian, has been suggested. The W.R.L.F. equations in polar- 
coordinates [2] relate power mismatch with voltage corrections as 


■ AP ■ 
AQ 



■ AS 
AV/V 


( 2 . 12 ) 


Uhen 0 
Jacobian [J] 
is the number 
violation of 
[ J ] of size 
slack) as P-Q 
flow easily 


-limits of the sources are considered the size of 
will be (2M-Mg + ffl -1) x (2N - + m - 1), where m 

of P-V buses converted to P-Q type following the 
generator Q~limits. Consider an extended Jacobian 
(2N-2)x(2N-2) considering all source buses (except 
type. The [J ] can be formed at the end of load 
by augmenting SQ/SS and CSQ/SV.V) elements 
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cot'C'ftspondiri^j to &11 P— V buses in the final Jacobi an [ J 3 • A 
aenaitivity matrix [S] can be defined as [S] = vhich is 

formed at the end of base load flow directly provides sensitivity 
relations between bus powers and voltages and can be used to 
compute new bus voltage angles and magnitudes, if the changes in 
bus power injections are known. The relationship can be written 
as 

■ A<5 ■ 

AV/V 


[^] 


AP 

AQ 


C2.13) 


If the generator or line outages are simulated as changes in 
bus real and reactive power injections, the post outage changes In 
voltage magnitude and angles can be computed using Eqn. (2.13). 
Uith the known new complex bus voltages, the new reactive power 
outputs of sources and hence the changes in the outputs can be 
computed. The attractive feature of this approach is that it does 
not require any additional load flow iterative simulation of 
contingencies and hence distribution factors and post outage 
conditions can be calculated very fast. The procedure for 
simulating line and generator outages in the sensitivity 
relationship and clauclation of the distribution factors is given 
below. 


2. A. 1 Li ne Outage Voltage Distrlbuticin Factors 

Fig. 2.2Ca) shows the pre-outage state of a part of a power 
system network, where line— i conencting bus i and j is to be 
considered for outage study. 
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Fig. 2.2(b) ahowa the post outage state of power systeia in 

which liiie**“l is out of service. Usual simulation of the line 

outage reQuires modification of exclude the parameters 

of line*“l which changes the Jacobinan and hence will involve a 

time extensive process. In order to retain the original *Y 

^ bus " 

elements of Jacobinan and sensitivity matrix, line outage has 
been simulated by considering two fictitious generators at bus-i 
and bua-“jy the two terminal ends of the line. If the power output 
of the two fictitious sources are same as the line flows at the 
two ends, the net power will be zero thus simulating the outage 
condi t ion . 

Thus, changes in bus powers from pre-outage to post-outage 
state at buses i and j for outage of line-l are, 



AP . 

= P. . 

. AQ, = Q. . 
i i j 





1 

i J 





AP . 

= P . . 

and AQ . = Q . . 



(2.14) 


J 

J 1 

J JX 



For 

the outage of 

the line-l these 

four elements 

are 

entered in 

the 

right 

hand side vector of Eqn. 

(2.13) and all 

other 

remaining 

el ement a 

in this 

vector will be zero. If either 

bus-i 

or J is the 


slack bus, only two non-zero elements will appear in the [AP , AQ] 

vector - 


The solution of Eqn. (2.13) gives the changes in bus voltage 
angles and magnitudes from pre-outage to post-outage condition. 


Uith the changes in 

voltage 

magnitude at all 

the 

bus es 

known 

(alack bus voltage 

assumed 

to be constant), 

the 

1 ine 

outage 



voltage distribution factors can be computed using Eqn . (2.4) 

T 

where is computed from base load flow results using Eqn. (2.3). 
The factors can be calculated for each line outage, considering 
one at a time for (I = l,2^..,Nj^), on similar approach. Thus the 

total number of LOVDF ’ a will be N«xN, which can be stored in 

s L 

computer memory, say aa matrix [A], to carry out voltage 
contingency analyaia. 


2* 4* 2 Gemerator Oixtage? Voltage Dlslrlbullon Factors 


•Fig, 2.3(a) ahowa the part of 
generator g, which La connected to hue i, 
for outage atudy. Fig. 2.3(b) ahowa the 
power ayatem in which generator g ia out 
a t at e 


a power ayatem, where 
ia under consideration 
post outage state of 
of service. In this 



and Oi = Oi - 0^ 

where P| and are respectively the real and reactive 

bus-i in post outage state. 


(2.15) 


power of 


Thus changes in ith bus power from pre-outage to post-outage 
are given by 

AP, = - P and AQ . = - (2.16) 

i g 1 g 

All other elements of AP and AQ vector will be zero. Uith these 
AP and AQ , the Eqn. (2.13) is solved to compute changes in bus 
voltage magnitudes (AV) . The GOVDFs can be computed for outage of 
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each generator taken one at a time using Eqn . (2.6). The GOVDFs 
thus obtained, are stored in the computer memory, say as matrix ■ B 
for conducting contingency analysis. Outage of slack bus 
generator can not be simulated in Eqn. (2.13) as it excludes 
elements corr eaponding to the alack bus. However, its simulation 
will require considering another generator bus as alack and 
redefining matrix [S] accordingly. In the present study, outage 
of the alack generator has not been considered. 

2.4.3 Li tre Outage Reactive Povrer Dlstribu-tlon Factors 

In certain voltage security studies, the effect of outages 
on outputs of reactive power sources, instead of bus voltages are 
required to be computed. The calculation of such factors requires 
the base load flow to be run without considering the Q-limits of 
sources. Thus the sensitivity relationship is redefined as, 


■ A<5 ' 

- r s’l 

■ AP ■ 

_ AV/V_ 

■ J 



where [S'] is the inverse of load flow Jacobian directly available 
at the end of base load flow and has dimension 
(2N-N^-1 )x(2N-Nq- 1) . Consider the outage of a line-l. The change 
in bus voltages are calculated exactly in the similar lines as in 
section 2.3 for voltage distribution factors but using Eqn. (2.17) 
instead of (2.13). The LOQDFs can be computed using Eqn. (2.5) 
for each line outages simulated one at a time. In this case if 
either bus-i or j is P-V type bus then the corresponding AQ 
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element will not appear on RHS vector in Eqn. (2.17) thus reducing 
the non-zero elements to 3. 

Using the computed new voltages and the bus reactive powers 
at Q-aources and hence their reactive power output and the change 
in reactive power outputs from pre-outage to post— outage can be 
computed. LOQDFs can be calculated using Eqn. (2.5) for each line 
outage taken one at a time. 

Thus, the total number of LOQDFs is Nj xN„ which are selared 

b W 

in the memory of computer, say as matrix [C] for contingency 

analysis. 

2. 4. 4 Gernorat.or OuLager Reactive Power Dlstrlbvrtlon Factor’s 

Consider the outage of a generator-g. The change in bus 
voltages magnitude and angles are computed exactly on the same 
lines in section 2.4.2 but using Eqn. (2.17). The post-outage 
voltages are thus computed by adding the changes in voltages to 
the base case voltages. Using the post-outage voltages, the bus 
reactive powers at all other source buses and hence the reactive 
outputs of the Q-aources and the change in reactive power 
generation from pre-outage to post-outage are computed. The 

GOQDFs can be calculated for outage of all the generators taken 
one at a time using Eqn. (2.7). These GOQDFs are stored in the 
memory of the computer, say as matrix [D] for contingency 

analysis. 
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2.5 SYSTEM STUDIES AND RESULTS 

To 6£it:abli8h the e£ f ect iveneas of the proposed distribution 
factors, studies were carried out on two sample systems : (i) 
IEEE-14 bus system as described in Appendix-A and (ii) a practical 
19-bus Indian system representing 400 kV, UP State Electricity 
network described in Appendix-B. The IEEE-14 bus system contains 
20 lines, three generators and two synchronous condensers whereas 
the Indian system contains 24 lines and 4 generators. 'For the 
present study, contingencies of all the lines taken one at a time 
for both the systems and outage of only two generators (gen-2 & 3) 
in 14 bus and outge of three generators (gen-2, 3 & 4) in 19 bus 
Indian system, taken one at a time, were simulated. 

The computer simulation of both the systems were carried out 
on HP 9000 computer. The distribution factors were computed for 
base case loadings as mentioned in Appendices A and B. The line 
generator outage voltage and reactive power distribution factors 
for IEEE-14 bus system are given in Table 2.1. To establish the 
accuracy of the distribution factors method, for obtaining 
post-outage bus voltage magnitudes and the reactive power outputs, 
the results have been compared with those obtained using full AC 
load flow method employing NRLF in polar coordinates. The 
comparison for some of the contingency cases are presented in 
Tables 2.2 and 2.3. It can be observed that the bus voltages 
obtained using the proposed distribution scheme are very close to 
that of full AC load flow with a error of around 3%. However, the 
error in predicting the reactive power output is more. 


The full 



AC load flow did not converge for line~l outage, whereas the 
results of post"“OUtage voltages were obtained through distribution 

factors., 

The voltage and reactive power distribution factors coxnputecL 
for the 19 bus Indian system are shown in Table 2.4. The 
comparison of post-outage bus voltages and reactive powers using 
full AC load flow and distribution factors technique for few 
contingency cases are shown in Tables 2.5 & 2.6. It is 

interesting to note that out of the 27 contingencies simulated, 
the full AC load flow method converged only for 15 cases and did 
not converge for remaining 12 cases. However, the post outage 
state can be predicted using the distribution factors method. The 
comparison of voltages (Table 2.6) shows that the voltages 
obtained from distribution factor methods are quite close to that 
obtained by exact load flow. The error is around 6%. The error 
in predicting the reactive power is however significant. 

Further, it was felt worth exploring the accuracy of 
predicting post outage states for changes in system loading. To 

study this effect, the loading of 14 bus system was increased by | 
51 and 10% and was decreased by the same 5% and 10% from their 
base values. The real and reactive power loadings were changed at | 
all the buaes aimultaneoualy to aimulate the above condition. The j 
compariaon o£ voltage predicted uaing the dlatribution factora and | 
the full AC load flow method for a aample contingency (outage of 
line~J4-) are preaented in Table 2.7. It can be observed from this 
table (Table 2.7) that for decrease in loading even upto 10%, the 
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prticlictfed voltage values are cjuite accurate. However, for 
increased loading condition the results are accurate only upto 5% 
increase. Thus, the distribution factors need not be recomputed 
for alight change in loading in this range. 

a. 6 CONCLUSIONS 

In this chapter new set of voltage and reactive power 
distribution factors for both line and generator outages have been 
suggested. The voltage contingency studies conducted on two 
sample systems reveal that, 

i) The prediction of post outage bus voltages using the 
distribution factors are quite accurate and provides the 
results with a error of 3% in IEEE-14 bus system and 5% in 
the Indian system. However, with the use of distribution 
factors, the error in predicting the reactive power output of 
sources is higher. 

ii) Since the distribution factors are obtained directly from 
base load flow results without involving any iterative scheme 
or the additional load flow solution, its calculation and 
updating is quite fast and can be carried out on real time 


basis . 
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iii) The set of dlatrlbution factors computed at a base loading 
accurately predicts the post outage voltages of the system 
even for small change in system loading. Thus, they need not 
be recomputed for small deviation in system loading (about ± 
5%) This further reduces the computational time for the 
voltage contingency analysis. 
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TABLE-a. la 


LINE OUTAGE VOLTAGE DISTRIBUTION FACTORS (LOVDFs) 




- 







M ^ .MM* 







< -- 

— 


— Bus No 

— 



y 

Line 

No 


1 


2 

3 

4 


5 

6 

7 

1 



.000 


031 

. 029 

. 031- 


. 030 

. 031 

. 032 

2 



.000 


S95 

1 . 340 

. 956 

1 

. 255 

1 . 273 

1.313 

3 



.0 00 


1 31. 

. 23 9 

.614 


. 322 

. 326 

. 335 

4 



000 


028 

-.511 

- . 222 


. 442 

- . 448 ' 

-.466 

5 



000 

-2 , 

320 

-17. 959 

-11. 528 

-20 

. 564 

-20 . 352 

-21 .402 

6 



0 00 


0 08 

. 070 

-.018 


. 080 

- . 0 82 

- .1 08 

7 



000 

- 

005 

-.013 

-.010 

- 

. 024 

- . 024 

-.030 

3 



00 0 


006 

. 004 

- . 012 

- 

. 037 

- . OIB 

- . 0 45 

3 



0 0 0 

- 

0 1 8 

“* , 033 

- 015 

- 

. 036 

- . 037 

- 0 37 

1 0 



0 00 


00 1 

.0 02 

.001 


.171 

. 004 

. 0 03 

1 1 



0 0 0 

- 

0 07 

- 04 1 

-.011 


. 001 

.00 1 

-.103 

1 2 



0 00 


000 

- 020 

.001 


.012 

.012 

.017 

1 3 



0 00 


OO 1 

- 072 

.000 


. 028 

. 028 

. 041 

1 4 



00 0 


0 03 

-.134 

. 002 

- 

. 056 

- . 056 

-.091 

1 S 



.000 


0 02 

- . 050 

,0 06 

- 

. 122 

-.124 

- . 081 

1 6 



.000 


0 1 2 

- . 060 

. 023 

-- 

.015 

-.015 

- . 059 

1 7 



0 00 


03a 

- . 034 

- . 053 

- 

. 055 

- . 086 

- . 089 

1 8 



000 


01 1 

-.101 

. 024 


. 1 06 

. 1 08 

. 1 40 

1 3 



0 00 

- 

002 

. 015 

- . OD'4- 

- 

. 01 3 

-.013 

. 0.25 

20 



.000 


01 7 

. I oa 

- . 036 


. 1 45 

-.147 

-.190 





/ — 

— 


— Sus No . 

— 



) 

Line 

No 


a 


9 

1 0 

1 1 


1 2 

1 3 

1 4 

1 



. 0S7 


028 

.032 

. 031 


. 030 

. 030 

. 032 

2 


1 

.260 

1 . 

1 75 

1 . 327 

1 . 335 

1 

. 361 

1.361 

1.361 

3 



. 256 


308 

. 331 

.313 


.237 

. 303 

.328 

4 



. 433 


404 

- . 476 

-.494 

- 

.517 

-.514 

- . 436 

5 


! 5 

, 683 

-13. 

745 

-21 . 043 

-13.531 

-18 

. 472 

-18.311 

-20 . 777 

6 



. 005 

• ^ 

026 

-.163 

-.316 


. 053 

. 041 

- . 046 

7 



003 

* ^ 

01 3 

- , 028 

- . 022 

- 

. 488 

-.120 

- . 070 

6 



. 009 


0 1 6 

- . 040 

- . 020 

- 

. 033 

-.235 

-.131 

9 



. 027 


036 

- . 036 

- . 035 

- 

. 034 

- . 034 

- . 036 

1 0 



,001 


002 

.0 03 

. 0 03 


. 002 

. 002 

. 003 

1 1 



... 0 1 a 

_ 

0 1 3 

- . 034 

- . 070 

- 

. 045 

- . 053 

- . 083 

1 2 


»- 

00 1 


002 

- . 070 

- . 046 

- 

.017 

-.015 

. 003 

1 3 



.010 


00 1 

023 

- , 022 


. 1 02 

-.133 

- . 338 

1 4 



003 


0 02 

- 1 03 

- 1 22 

- 

1 32 

- 134 

_ 117 

1 5 



0 03 


■1 0 P 

^ 97 J. 

-.061 


.05*^ 

- 05 4 

-.070 

1 6 



0 2: 


: 3 = 

- ^ 

- . 052 

- 

0 6 4 

— . 0 5 0 


1 7 


«*» 

OcO 


O -7 

0 0 

. w w 

-.061 

- 

. 037 

- , 0 44 

^ 72' 

1 3 



0 02 


J S "Z 

227 

- . 333 

- 

. 053 

-.06 0 

:} 3 7 

1 3 



: 0 0 2 

- 

006 

- : 0 1 9 

- . 0 03 
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-.147 

- . 08 0 

SO 



013 

* 
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- 148 

035 

— . 

.171 

.231 
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TABLE-E, lb 

GENERATOR OUTAGE VOLTAGE DISTRIBUTION FACTORS (GOVDFs) 


C#n . 

No . 1 

\ — — 

2 

3 

' Bus No - 
A 

5 

> 

6 

7 

E 

3 

,000 
. 000 

- . 069 
-.075 

-.061 
- . 362 

-.066 
- . 121 

- , 064 
-.211 

-. 064 
-.214 

- . 066 
- . 241 




/ — 






Gen . 

No , 8 

9 

1 0 

“ Bus No * 
1 1 

12 

13 

1 4 

2 

3 

- . 055 
-.156 

-.061 

-.157 

-. 066 
- . 263 

- . 064 

-.311 

- . 063 
-.359 

- . 063 

- . 350 

- . 066 
- . 294 


TABLE &.1-C 


LINE OUTAGE REACTIVE POWER DISTRIBUTION FACTORS 

(LOQDFs) 


Line No . 



- Source 

3 




1 

2 

0 . 

4 

5 

1 

. 3770 

-1 .1835 

- . 0001 

- 

. 0006 

. 0000 

2 

-3.1415 

-.2444 

. 0004 


.0190 

,0004 

3 . 

-.9942 

.6615 

-.0006 

-1 

.3051 

. 0001 

4 

2.7673 

. 9050 

- . 0004 


. 0206 

-.0002 

S 

-8.7072 

. 0086 

-.0131 

- 

. 1601 

. 001 1 

6 

. 1901 

. 0000 

1 . 0926 


. 0000 

. 0000 

7 

. 1257 

. 0000 

1 . 0603 


. 0000 

. 0000 

0 

. 1395 

. 0000 

1 . 0693 


. 0000 

,0000 

9 

.2307 

. 0064 

. 0000 


. 9960 

. 0000 

1 0 

-.0104 

. 0000 

. 0000 


. 0000 

-1.0135 

1 1 

.•1817 

. 0000 

,0020 


.0000 

.0000 

12 

. OOtT 






1 3 

. 0551 

.0000 

.0012 


. 0000 

.0000 

1 4 

-.1177 

0002 

-1 . 1 995 


. 0000 

.0000 

1 5 

- , 0905 

. 0000 

-.0018 


. 0000 

. 0001 

1 6 

- . 5397 

- . 0001 

- . 0095 


.0000 

. 0000 

1 7 

. 2702 

. 0069 

.0016 


.0016 

.0000 

1 8 

- . 8045 

- . 0002 

- . 0307 


. 0000 

. 0000 

1 9 

0434 

.0000 

. 0025 


. 0000 

. 0000 

20 

.5172 

. 000 1 

. 0145 


. 0000 

. 0000 


TABLE 2. Id 

GENERATOR OUTAGE REACTIVE POWER DISTRIBUTION FACTORS 

(COQDFs) 




- Source No. 

3 



< 

C»n. No. 1 

2 

4 

5 

2 1.1 904 

. 0000 

.0004 

. 002S 

-.0001 

3 1 .8875 

. 01 08 

.0000 

.0009 

-.0004 



TABLE-a.£ 


COMPARISON OF POST-OUTACE VOLTAGES OBTAINED USING FULL AC LOAD 

FLOW AND DISTRIBUTION FACTORS METHOD 

Bus No , 

L i ne- 1 0 

ACLF Factors 

Sen 

ACLF 

- 2 

Factors 

Line~15 

ACLF 

( Ira ns former- 8 ) 
Factors 

1 

1 , 000 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1.000 

a 

0 . 990 

0 , 992 

0.961 

0.958 

0.991 

0 . 998 

3 

1 .000 

1.000 

0 . 991 

1 . 000 

1 . 000 

1.006 

4 

0 . 965 

0 . 969 

0 . 943 

0 . 969 

0 . 967 

0 . 968 ► 

S 

1 000 

1.014 

1 . 000 

1.000 

1 . 000 

1.014 

6 

0 . 972 

0 . 986 

0 . 975 

0 . 986 

0 , 983 

1 . 000 

7 

0 . 970 

0 . 979 

0 . 967 

0.980 

0 . 972 

0.989 

8 

0 . 964 

0 . 967 

0 . 949 

0 . 967 

0 . 965 

0 . 967 

9 

0 . 958 

0 . 962 

0 . 942 

0 . 962 

0 . 960 

0.961 

1 D 

0 79157 

0 . • 9 TS 


4 . 97S 



1 1 

0 . 979 

0 . 983 

0.973 

0 . 983 

0 . 979 

0 . 990 

1 2 

0 983 

0 . 984 

0 . 974 

0 . 984 

0 . 983 

0 . 990 

t 3 

0 . 976 

0 . 978 

0 . 968 

0 . 978 

0 . 976 

0 . 984 

1 4 

0 , 953 

0 . 959 

0.948 

0 . 959 

0 . 954 

0 . 967 


TABLE-2.3 


COHFARISON OF 

POST-OUTAGE Q 

-SOURCE 

REACTIVE 

POWER OUTPUTS 

OBTAINED 


BY FULL 

AC LOAD FLOU 

AND DISTRIBUTION 

FACTORS METHOD 

I-Source 

Line-10 

Gen . - 

'2 

Line-15 ( Transf ormer-i 

No. 

ACLF 

Fact ors 

ACLF 

Fact ors 

ACLF 

Fact ors 

1 

-0 . 4480 

-0.4524 

-0 . 0994 

0.1539 

. 4477 

— t/ . 444 4- 

2 

0,5867 

0.5098 

0 . 0000 

0 . 0000 

0.5140 

0.5092 

3 

0.1673 

0.1440 

0.1980 

0 . 1 442 

0. 1462 

0.1441 

4 

0 . 6491 

0 . 6392 

0 . 8231 

0 . 6371 

0 . 6378 

0 . 6382 

5 

0 01 9£ 

0.1017 

0.0771 

0 . 0520 

0 . 0873 

0 . 0521 



TABLe-2.4a 


LINE OUTAGE VOLTAGE DISTRIBUTION FACTORS (LOVDFs) 


nt No . 

1 

1 

— 





2 

3 

— - DUS NO . 

4 

5 

6 

7 

1 

000 

055 

. 065 

. 059 

. 063 

. 051 

065 

2 

000 

• . 045 

- . 046 

-.038 

-.044 

027 

- . 046 

3 

0 0 0 

. 029 

036 

-.027 

- . 029 

023 

- . 027 

4 

000 

- . 02^ 

- . 024 

-.035 

-.023 

-.017 

- . 024 

5 

, 000 

- . 00^ 

- . 005 

- . 004 

-.004 

-.003 

005 

6 

.000 

-.743 

-.738 

-.676 

-.736 

-.561 

- . 741 

7 

000 

62^ 

. 629 

.569 

.619 

.472 

.631 

8 

000 

.011 

.012 

.012 

.01 1 

.010 

.012 

9 

. 000 

. 026 

. 026 

.023 

. 026 

. 018 

. 026 

t 0 

0 00 

.209 

.212 

. 183 

.207 

. 1 49 

.213 

1 1 

000 

.13? 

. 138 

.119 

.134 

. 096 

. 138 

1 Z 

000 

36? 

. 371 

.323 

.362 

.265 

. 372 

13 

000 

- ,28$ 

- . 289 

-.250 

- .283 

- .21 3 

- . 290 

1 4 

000 

-12. 08| 

-12.244, 

-10.813 

-11.963 

-8.079 

-12.291 

15 

000 

-.230 

- . 232 

-.208 

- .228 

-.171 

-.833 

1 6 

000 

1 

. 143 

. 128 

.141 

. 1 05 

144 

t 7 

000 

26f 

271 

.242 

.266 

. 1 99 

. 272 

1 e 

000 

294 

, 296 

264 

.291 

.217 

. 297 

1 9 

0 00 

51? 

-.517 

-.460 

- .507 

- . 378 

-.519 

20 

000 

S85 

- , 508 

-.528 

- .579 

- . 429 

- 590 

21 

000 

- 0 1 9 

-.019 

-.017 

-.010 

-.013 

- 01 9 

22 

000 

• t S3 

-. 153 

-.138 

-.151 

-.114 

- . 154 

23 

000 

-1 31$ 

-1 320 

-1.194 

-1 .303 

- . 985 

-1 . 325 

24 

000 

- . 252 

- .254 

-.227 

- .250 

187 

- 255 


( ^ « 


Lin# No 

8 

9 

1 0 

1 

060 

08? 

. 089 

2 

•* 0 38 

■■ OSS 

- 055 

3 

* ,087 

- 0 37 

~ . 0^0 

4 

- 080 

- 0 88 

031 

5 

« 0 04 

- 005 

- . 005 

6 

- 680 

- 880 

- . 936 

7 

,578 

740 

.787 

8 

.018 

,017 

.018 

9 

083 

030 

. 035 

1 0 

, 1 84 

140 

.808 

1 1 

. 180 

. 096 

.113 

12 

. 385 

381 

,384 

13 

- . 859 

- . 894 

-.285 

1 4 

10,869 

-18.843 

-12.536 

IS 

-.809 

- .846 

-.849 

16 

189 

130 

. 158 

17 

843 

.885 

.890 

18 

.265 

.3^9 

.313 

19 

» 463 

-.539 

-.551 

80 

- 531 

- 678 

-.717 

81 

- , 0 1 7 

- oeo 

-.081 

22 

- , 1 39 

- 178 

-.189 

83 

-1 .200 

-1 539 

-1 ,688 

84 

»- , 288 

-.869 

- .274 

... ,, 

... ... „ 


- -- - — 


Bus No . > 


1 1 

12 

13 

1 4 

.118 

.163 

. 1 68 

, 1 65 

-.072 

-.105 

-.108 

'.107 

050 

-.073 

- . 075 

- . 074 

- . 040 

-.050 

- . 060 

- , 059 

-. 007 

-.010 

-.010 

-.010 

-1.148 

-1 .587 

-1 . 620 

-1601 

. 961 

1 . 340 

1 . 376 

I 354 

. 081 

. 031 

. 031 

. 031 

. 043 

. 062 

. 064 

. 063 

. 887 

.443 

. 456 

. 451 

.178 

.874 

.283 

. 28 0 

,478 

.746 

.769 

. 759 

-.317 

- .584 

-.541 

- 536 

-14.842 

-83.670 

-24.417 

-24 . 147 

-.867 

-.347 

-.361 

- 360 

. 166 

155 

. 166 

174 

.321 

.380 

.301 

.303 

. 343 

.386 

.345 

. 352 

-.627 

-.763 

-.777 

- . 757 

-.879 

-1 .223 

-1 .856 

-1. .837 

025 

-.033 

- . 034 

- . 033 

- . 228 

-.301 

-.307 

- . 300 

-1 . 966 

-2.597 

-8.647 

-2 . 583 

-.305 

-.308 

-.303 

- 283 






cont d 



Line No . 


16 17 10 


1 .114 ,114 

£ . 075 - . 07S 


3 

- . 061 

050 

4 

- . 041 

041 

5 

- . 007 

007 

6 

-1 . 043 

-1.034 

7 

.886 

.879 

8 

, Ogl 

. 021 

9 

. 044 

. 043 

1 0 

. 334 

.333 

1 1 

.214 

,213 

12 

.571 

.570 

1 3 

- . 426 

-.426 

14 

-18. 569 

-18.535 

15 

-.318 

-.318 

1 6 

. 1 95 

. 1 95 

1 7 

,364 

.366 

18 

. 402 

.404 

1 9 

- . 688 

-.692 

EO 

- 807 

- . 806 

ai 

- . 0^0 

- . 024 

22 

- . 182 

-.100 

23 

-1 , 568 

-1 . 547 

24 

- . 340 

-.342 


.143 

.145 

. 165 

-.094 

-.095 

-.107 

-.063 

-.064 

-.074 

-.051 

-.052 

-.059 

-. 009 

- .009 

-.010 

-1 .297 

-1.315 

“1 .600 

1.102 

1.118 

1 .353 

. 026 

. 027 

. 031 

. 054 

. 055 

. 0'63 

.418 

.424 

.450 

.268 

.271 

.280 

.714 

.724 

.759 

-.534 

-.542 

-.535 

-23.237 

-23.564 

-24. 130 

-.399 

-.404 

-.360 

.245 

.248 

. 174 

.458 

.465 

.-303 

.506 

.513 

.352 

-.867 

-.879 

- .757 

-1.010 

-1 . 024 

-1 . 236 

-.030 

-.031 

- . 033 

-.177 

-.179 

-.300 

-.553 

-.482 

-2.581 

-.428 

-.434 

-.305 


TABLE-2. 4b 


GENERATOR OUTAGE VOLTAGE DISTRIBUTION FACTORS (GOVDFs) 




< 

1 


Bus No . 
3 



Qen , 

No 

a 

4 

5 

2 

3 

4 


. 0096*^00 
. 00(jE+00 

oop+oo 

-.146E*»-07 - 

, 1 48E-0t 
. I23E-01 

. 170E-f07 
.788E-0a 
. 126E-01 

-.160E+07 - 

. 133E-0t 
-.630E-03 

1 66E+07 
. 147E-01 
. 1E2E-01 










<• 

6 


Bus No. 
8 


\ 

Gtn . 

No 

7 

9 

— / 

1 0 

2 

3 

4 

— 

. 1 SQE + 07 

1 OlE-01 
, 781 E-0£ 

-.171E+07 - 

. 141E-01 
. 126E-01 

. 16tE+07 
. 134E-01 
. t 02E-01 

-.229E+07 ^ - 
. 175E-01 
. 1 44E~0 1- 

. 257E + 07 
. 180E-01 
. 1 49E-01 












Bus No . 
13 


\ 

Cen . 

No 

< - 

1 1 

ia 

1.^ 

/ 

15 

a 

3 

4 

- , 

.30SE+07 

.244E-01 

ao9E-oi 

-.439E+07 - 

.361E-01 
.295E-01 

.45IE+07 

.373E-01 

.303E-01 

-.443E + 07 
.369E-0r 
.299E-01 

. 304E + 07 
aeoE-oi 
,21 OE-01 







>■ , 




Du. 

Klrt _ — — . 

\ 


G#n . 

No . 

i - 

1 6 

17 

riO . 

18 

/ 

1 9 


2 

3 

4 

— . 

30()E+07 
869E-01 
£1 QE-01 

-.377E+07 

.326E-01 

.263E-01 

382E+07 

331E-01 

267E-01 

-,443E+07 
. 369E-01 
.299E-01 



TABLE a. 4c 


LINE OUTAGE REACTIVE POWER DISTRIBUTION FACTORS 

(LOQDFs) 


t X n e No 




1 

u oource No . 

2 3 

> 

4 

1 

9 . 8272 

-5.5705 

.0005 

. 0008 

a 

a . 71S7 

.0002 

.0004 

. 0003 

3 

t 1 . 8277 

, 0002 

-6 . 0433 

.0008 

4 

6 . 071 S 

. 0002 

.0 0 03 

-4 . 1 999 

5 

1.1437 

. 0000 

.0 00 0 

.0000 

6 

32.3986 

-.0589 

-.0285 

-.0070 

7 

-30 . 1851 

.2818 

. 1720 

. 0634 

a 

. 1 053 

. 0000 

.00 00 

.0000 

9 

2 5954 

. 0001 

.000 1 

. 000 1 

1 0 

-16, 2504 

. 0043 

. 00 06 

- . 0013 

1 1 

-6 . 8581 

-.0001 

- . 0005 

- .0005 

1 a 

-46 8585 

. 0213 

. 0063 

- . 0028 

1 3 

1 0 . 7590 

-.0015 

. 0002 

. 0009 

1 4 - 

154.6911 

.3120 

. 1 648 

.0367 

1 5 

1 1 . 7321 

-.0262 

-.0138 

- . 0034 

1 6 

-8 . 2976 

.0016 

. 0000 

- . 0007 

1 7 

- 1 5 9 1 1 7 

. 0047 

. 0008 

- . 0012 

1 8 

-16. 6533 

.0013 

0008 

- . 0015 

1 9 

24,69 1 8 

-.0831 

- . 0467 

-.0135 

ao 

22.8610 

-.0336 

-.0162 

- . 0029 

ai 

1 . 4885 

. 0000 

. 0000 

. 0000 

£2 

8 . 4967 

0093 

0041 

-. 0005 

23 

56 . 1 069 

-.0386 

-.0144 

.0010 

24 

14,2014 

. 0004 

. 000^ 

- 0005 


GENERATOR OUTAGE 

TABLE £.4d 

REACTIVE POWER DISTRIBUTION 
(GOQDFs) 

FACTORS 



^ ^ 





— ) 

Gen No . 

1 

2 3 

4 

2 

341 9E+01 

.0000 -.4363E-04 

. 3827E-04 

3 

.2291E+01 

-.1098E-04 .0000 

. 1 1 99E-04 

4 

161 3E4.01 

-,5923E-05 ~.6996E-05 

, 0000 


4 





TABLE- 

2.5 


COMPARISON OF 1 

POST-OUTACE 

VOLTAGES 

OBTAINED USING 

FULL 

AC 

LOAD FLOW AND DISTRIBUTION 

FACTORS METHOD 

•*’ *** 

* "** •* '**’ 





Su* 

No . 

Lln*-5 

G*n.-3 



ACLF 

Factors 

ACLF 

Factors 

1 


1 . 030 

1 . 030 

1 . 030 

1.030 

e 


1 . 030 

1 .015 

1 . 030 

1 . 034 

3 


1 . 030 

1 . 035 

1 . 072 

1 . 042 

4 


1 . 030 

1 . 033 

1 .030 

1 . 050 

5 


1 . 009 

1 .017 

1 .015 

1 .036 

6 


1 . 003 

1.016 

1 .007 

1 . 029 

7 


1 . OEO 

1 . oei 

1.012 

1 . 037 

a 


t . 007 

1.017 

1.013 

1 . 034 

9 


0 . 9S8 

0.985 

0.976 

1 . 007 

io 


0 . 9S4 

0.987 

0.975 

1 . 009 

11 


0.936 

0 . 984 

0.973 

1.014 

i2 


0 . 930 

1.014 

0 .998 

1 . 059 

13 


0 . 95E 

1 . 028 

1.012 

1 . 074 

14 


0.925 

1.019 

1 . 003 

1 . 065 

Is 


0.958 

1.010 

1 . 000 

1 . 042 

16 


0 . 960 

1.011 

1 . 001 

1 . 043 

17 


0 . 976 

1 . 042 

1 . 029 

1 . 082 

la 


0 . 970 

1 . 044 

1 . 032 

1 . 085 

19 


0 . 947 

1 . 021 

1 . 005 

1 . 066 


TABLE-E.6 

COMPARISON OF POST-OUTACE Q-SOURCE REACTIVE POWER OUTPUTS OBTAINED 


BY 

FULL AC LOAD 

FLOW AND 

DISTRIBUTION FACTORS 

METHOD 

Q-Sourc# 

Lina 


Can . - 

2 

No , 

ACLF 

Factors 

ACLF 

Factors 

! 

0.613 

0.347 

1 .254 

2.002 

2 

1 . 686 

1 . 034 

1 .281 

1 . 034 

3 

0.954 

0.645 

0 . 744 

0.455 

4 

0.627 

0.455 

0. 000 

0 . 000 
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TABLE-8.7 


COMPARISON OF VOLTAGES PREDICTED BY THE DISTRIBUTION FACTORS METHOD 
AND THE FULL AC LOAD FLOW METHOD FOR A SAMPLE CONTINGENCY 
WITH DIFFERENT LOADING CONDITIONS 


Bui» No 

0 . 9 
ACLF 

PU 

Factor* 

0.95 PU 

ACLF Factor* 

1 . 05 
ACLF 

PU 

Fact or* 

1.10 PU 
ACLF Facto 

1 

1 . 000 

1 . 000 

t . 000 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1 . OO 

E 

0 . 999 

0 . 999 

0.992 

0 . 996 

0.998 

0 . 998 

0 . 989 

0 . 98 

3 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

1 . 000 

0.993 

0.99 

4 

0 . 98E 

0 . 982 

0 . 969 

0.976 

0.961 

0 . 961 

0.952 

0 . 95 

5 

1 . 000 

1 . 000 

1 . 000 

1 .000 

1 .000 

1 . 000 

1 . 000 

1.00 

6 

0.992 

0.992 

0.986 

0.989 

0.983 

0.983 

0 . 977 

0 . 97 

7 

0 . 987 

0 . 987 

0.979 

0.983 

0.975 

0 . 975 

0.968 

0 . 96 

8 

0 . 975 

0 . 975 

0.967 

0.971 

0.963 

0 . 963 

0 . 957 

0 . 95 

9 

0 971 

0 . 971 

0.962 

0 . 967 

0.957 

0 . 957 

0 . 950 

0 . 95 

) 0 

0 . 982 

0 . 982 

0.975 

0 . 978 

0 , 971 

0 . 971 

0 . 963 

0 . 96 

! 1 

0 987 

0 . 987 

0 . 983 

0 . 985 

0.981 

0 . 981 

0 . 973 

0 . 97 

1 E 

0 969 

0 . 987 

0 . 988 

0 . 986 

0 . 988 

0 . 985 

0 . 980 

0 . 97 

! 3 

0 . 979 

0 . 981 

0 . 976 

0 . 979 

0 . 979 

0 . 975 

0 . 966 

0.96 

1 4 

0 . 966 

0 . 967 

0.958 

0 . 963 

0 . 959 

0 . 955 

0 . 995 

0 . 99 



CHAPTER - 3 


PERFORMANCE INDICES FOR VOLTAGE CONTINGENCY SELECTION 

3.1 I NTRODUCTI ON 

Cont 1 ng ei'iC y aelection problem la coricerned with developing 
computer algorithma for quickly identifying those contingencies 
which may cause out-of-limit conditions so as to reduce the number 
of contingencies that need to be evaluated by full A.C. load flow 
when assessing the power system's security. Almost all 
contingency ranking algorithms employ a scalar performance index 
(PI) which is a function of chosen network variables. For 
contingency analysis, the PI value represents the degree to which 
the system operating state is approaching or exceeding its static 
limits. 

In planning studies, the traditional approach for steady 
state contingency analysis is to test all contingencies 
sequentially to evaluate system performance and reliability. This 
analysis consists of simulating outages of one or more generating 
and transmission units to investigate their effects on bus 
voltages and line power flows. However, exhaustive contingency 
testing including all conceivable contingencies, becomes a 
prohibitively time extensive when all primary contingencies, each 
with additional levels of secondary contingencies are involved. 
On the other hand, applying contingency testing to a subset of 
contingency cases selected only on the basis of the planner a 
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experietice and intuition may be inadequate due to the possibility 
of omitting some critical cases. 

In real time, as the system conditions change, the 
contingencies which cause insecure operation may also change, and 
would be different from those predicted by off-line simulation 
studies. Therefore, the selection of contingencies to be studied 
by the on-line security analysis program should be based on the 
current operating conditions, but not a fixed list based upon 
off-line studies, Threfore, ranking based method for contingency 
selection is popularly used. The method involves calculation of 
post outage state of system for each contingency using fast but 
approximate methods and computing their relative severity and 
ranking theai using some performance index. 

For voltage contingency selection, various performance 
indices with bus voltages or reactive power outputs of sources as 
the monitored variables, have been suggested- They have 
invariably utilized the second order performance index, which, In 
general, suffer from 'masking* and 'misranking' effects. Some of 
the efforts in reducing these drawbacks include the works of 
Halpin [16] and Schafer [25] who have suggested the method to 
compute the threshold value and optimal weights of the second 
order performance indices. The method suggested in ref. [16] for 
optimal weight selection involves a optimization technique and is 
qui t e Involved . 

In this chapter various existing performance indices for 
voltage contingencies have been critically reviewed. A more 



44 


efective set of voltage and Reactive power indices have been 
escploted which are of higher orders, A simple and fast approach 
for optimal selection of weights have been suggested which 
alongwith the proposed higher order performance indices elimina.tes 
misranking and masking effects. The outage studies have been 
simulated using the distribution factors suggested in Chapter-2 
and the ef f ect i veness of the new performance indices have been 
demonstrated on IEEE-14 bus and 19 bus U.P, State Electricity 
network , 


3.3 EXISTING PERFORMANCE INDICES 


Use of performance index for ranking of cont ingencies 
according to their relative severities was first suggested in 1979 
by Ejebe & Uollenberg [6] for both line and voltage securities. 
Various modified versions for voltage contingency selection were 
suggested by Albuyeh et al . [10], Medicherla et al . [11] tand 
Uasley et al . [14], The nature of performance indices alongwith 
the choice of weighing factors are summarized in Tables JT. 1 and 
J^.2 [16] considering respectively the impacts of bus voltage 
magnitudes and injected reactive powers. 


The general 

^ = I “i 

i 


form of the various performance indicies 



is 


(3.1) 


where f (Z) is a linear function of Z. which denotes either 

i - . 

[aV . /AV or rAQG^/AQGV"'''^*j , the changes in load bus voltage 



magnitudes or generator bus injections with respect to their 
target value or ratings respectively. It can be seen from tables 


Table 1 s Summary of PI used f'or measurlniQI bhe Impact of 

contntirigericies ori the bus voltage magnitudes. 



Normal ized 
Variable 


Neighing 

Factor 



z , 

1 


w . 

1 



E j ebe , 

Uollenberg [6] 

AV. 

1 

AV^ iro 
i 

z . 

1 

€jL . 

1 






4^ . AV . 

1 1 

v» • 

if z ^ 

< -1 

Albuy eh , 

Bos e , 

Heath (101 

AV. 

1 

. im 

AV , 

1 


0 

1 

< z, < 

1 

1 




oi . AV . 

i i 

M ’ 

v? 

i 

if z ^ 

> 1 

Medi cher la , 
Rastogi [11] 

AV. 

1 


2P .AV^^“ 

1 1 


^vliro 

v° 




i i 


AV. 



z 


i 


0 , 


if 



1 


Uas 1 ey , 

Daneshdoost [14] 


1, if |2.1 > 1 
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Table? 2 


StiiBmary of PI used For Measuring the impact of 
contingencies on the injected reactive power. 



Woraialized 

Variable 


Ueighting 

Factor 


z . 

i 

f .(2) 

w 

i 


E j e b fc , 

Uo 1 1 enb et'g [ 6 J 

QG^/QG^'^“ 

Z . 

1 


i 




oi . AV . 
i 1 

V” • 

if z < -1 

Albuyeh , 

Bose , 

Heath [10] 

aqg./aqgJ^“ 


0 

1 

z . < 1 

1 




^ Im 
ca.AQG, 

i i 

H ' 

if z^ > 1 


Uas 1 ey , . . 

Daneahdooat [14] AQGj./AQG:^“ 



z . 
1 


1, if jz^l > 1 
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1 iind 2 that the major difference between the various Pis is 
the choice of performance function f(£) and weighing coefficient 
w^ . Some of the current works have also concentrated on the 
automatic contingency selection approach. The basic idea has been 
to compute the changes in PI (i.e., J) for each contingency using 
an efficient sensitivity method [24]. The work has concentrated 
on improving the accuracy of the prediction of J obtained via 
these first order sensitivity methods. 


Another work by K. Nara et al . [17] presented a new concept 
in formulating a PI for contingency selection concerning voltage 
security analysis. The performance index is defined as a second 
order vector norm in the voltage apace. Two types of voltage 
limits-alarm limit and security limit are introduced. Each axis 
of the voltage apace is normalized by the difference between these 
two limits. Contingency cases in which the performance indices 
are greater than one are selected as severe cases and can be 
ranked in order of severity. The performance index which Includes 
the effects of generator violations and voltage limit 
violations is defined as follows : 


PI 



d“CL) 12n 

e-CL) J 



. fi-CL) 


2n 


l/2n 


(3.2) 


where 



i 

From the PI 


the total number of buses in the system, 
the system condition is in, 
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a s e c u r e 

state 

if 

: PI 

= 0 

an 

s t at 

i f 

V 

o 

PI 


an insecure state if : PI > 1 

Though different PI ranking algorithais have been reported in 
the literature, but the common disadvantage of all these methods 
is, that computationally feasible evaluations have successfully 
been derived for an exponent of two only. Second order indices 
have a very limited ability to detect localized system overloading 
or voltage violations. As a result, severe localized overload 
lines or voltage violation cases can be given far less importance 
than would be judged necessary from an engineering aspect. This 
problem has been term^ed as the 'masking effect’ in contingency 
selection. The phenomenon of selecting soaie of the critical 
cont ingenci es as non-critical is defined as 'masking'. Another 
problem, which is mainly due to coaiputation of post outage state 
of the system using approximate methods during contingency 
slection is the 'misranking' which means inaccurate assignment of 
order of severity to various contingencies as should have been 
using the exact method. 

Recently, K.F. Schafer et al . have presented a new method to 
compensate the masking effect in second order PI algorithms [25]. 
For the PI distribution extended vector norm formulation is 
introduced, which additionally allows the quantitative description 
of the masking effect range. They defined the PI as follows 




49 


N 1 / 

2 (3.3) 

i = l 

where Is b . linear function of the system variables, which 

denotes e.g. (P. -/P- - or (AV./iiV. ), N the number of 

1 j 1 j , max i 1 , max 

regarded components (the number of buses and/or branches), m 
the exponent o£ the PI function, and w^^ a weighing coefficient. 
The subscript k should indicate the considered contingency k. 

In 1984, T.F. Halpin et al . [16] have presented a method for 

evaluating the effectiveness of the Automatic Contingency 
Selection (ACS) algorithms in capturing contingencies which give 
out-of-limit conditions. The evaluation is given in terms of both 
capture and flase alarm rates, which are computed apriori together 
with a stopping criterion. The Authors also pr.oposed a method 
which enhances the ACS algorithm’s effectiveness in terms of 
maximising the capture rate (CR) and minimizing the false alarm 
rate (FR). Theory has been given for selecting the set of weights 
in the scalar performance index for both real power and voltage 
problems in order to circumvent some of the misranking problems, 
and wore importantly, to select the threshold value of the PI in 
ordet- to capture all the critical contingencies while minimizing 
the capture of the non-critical ones. The capture rate and false 
alarm rate have been formulated using the probabilistic approach. 
The method is, in general, quite cumbersome and time extensive. 
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3 . 3 PROPOSaED PERFORMANCE INDICES FOR VOLTAGE CONTINGENCY RANKING 

Two aeta of performance indices are proposed to measure the 
relative severities of contingencies in terms of their effect on 
voltage levels and reactive power output of the sources. The 
volt^e level performance index chosen to quantify system 
deficiency due to out~of~lifflit bus voltages is defined as 


PI 


k 

V 


W 




vi 

n 


V^^n 


( ^^Hm ] 


( 3 . 4 ) 


where 


V 


sp 

i 


AV 


iim 


n 

N 


k 


voltage magnitude at bus 1 

specified (rated) voltage magnitude at bus i 

voltage deviation limit, above which voltage 
deviations are unacceptable 
exponent of penalty function 
number of buses in the system 
real weighing factor 
outage case number. 


The voltage deviation limit represents the threshold 
above which voltage level deviations are outside their limit. Any 
contingency case with voltage levels outside this yields a high 
value of the index PI^. On the other hand, when all the voltage 
level deviations from the rated voltage are within AV^ the 
voltas^ parfo«an=a indax la aaiall. Thus, thia index meaauree 
the severity of the out-of-limit bus voltages, and for a set of 
contingencies, this index provides a direct means of comparing the 
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relative severity of the different outages on the systetn voltage 

p V o file. 


It is pertinent to note, that since the bus voltage levels 
depend mainly on the reactive power flows, and therefore, on the 
reactive power production of the generators and reactive power 
sources such as synchronous condensers etc., the performance index 
PI^ provides a good measure of the severity of abnormal voltages, 
as long as the generating units remain within their reactive power 
liftiits.. However, it is possible to encounter a contingency for 
which some generator reactive powers are driven beyond their 
operating limits. Therefore, in order to reflect the reactive 
power capability constraints of the generators in the contingency 
selection for voltage analysis, we define another performance 
index Pip as follows 


OG.- Qg“''«^‘ 

1 X 


( 3 . 5 ) 



number of reactive power sources 
exponent of penalty function 
reactive power output of ith source 


= I [a=i,„ax - 0=1, «in] 
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III the above expressions QG , 

i,max ^ i,inin 


are 


and QG 

respectively the maximum and minimum limits of QG , The reactive 

1 ‘ 

power deviation limit AQG^^“ represents the span that is allowed 

for QGj^to move away from the target point in either direction 

(towards QG, or QG , , ) 

^ ^ i,atax ^ i,min'^' 


Contingency with the reactive power outputs outside this 

limit results in a greater, value of the other hand, when 

QGs are within the deviation limit AQg!^“ then the PI^ is small. 

1 Q 

IhuSythis index measures the severity of Q-limits violation of 

generating units for a particular contingency. 

The exponent ''n’' in the definition of the two performance 
indices and Pig in Eqns . (3.4) and (3.5) have been varied from 

2 to 20 to explore its desired value or the order of the 
performance indices , which minimizes the 'Masking effect*. For 
this set of study all the weights w^ were assumed to be one. In 

order to eliminate the misranking effects, the proper selection of 
Wj, is required. A simple and effective scheme is described in the 
subsequent s ec t i on . 


3,4 PROPOSED SCHEME FOR WEIOHTS ADJUSTMENT 

An important aspect in the use of P.I. based ranking methods 
is the proper selection of weights. The improper selection of 
weights causes misranking of contingencies which is mainly due to 
the inaccuracies of approximate load flow or distribution factor 
models in computing post outage states of the system. Some 
approximate guidelines for selecting the weights are found in 
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literature as given in Tables 1 and 2. However, as proper 
attempt has only been made by Haipin et al . [16] who used an 

opt Imi siat I <>«> technicjue t© compute the weights. The method 
uiaximizea the 'capture rate’ and minimizes 'false alarm rate' both 
formulated as pobabilistic functions and is quite involved. 
Hence, a simple approach has been proposed in this section, which 
is based on 'Least Square Error’ principle. 

The error is defined as the performance index value computed 

by the distribution factor (D.F.) method and the exact A. C. load 

flow method. If for kth contingency the performance indices 

(either Ply or Pig) from the two methods are fres p e c ti v e 1 -y 
k k 

(PI )0 p (PI ) gjj respectively the error is defined as 


E. ^fpi^^l - fpi^l 
^ -^DF ^ ^EX 


k € all contingencies 


(3.6) 


The proposed method minimizes a function F defined as the 
sum of the squre of errors 

F - J (E^)^ <3. 7) 

k = l 

where is the total number of contingencies to be considered. 

The optimal values of weights are computed off line at some 
average base operating points and can be used for subsequent 
calculation. The performance indices computed at this base 
operating point for different contingencies using exact lecxd flow 
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frt o d f:J 1 c an b e 

ti s r d, €: e i d e cl 
k 

probl eiii , ( P I 
to be^ equal 


coiifiputed cons iderifi^ all weight ages as unity or the 
values, if any. Thus, for the above minimization 
^ EX apriori known and constant. Let us assume it 

to (k =1, . . 


The post outage states and hence the value of V, AV, . . , QG 

1 imi t 

and j: .|- Eqns . (3.4) and (3.5) will be known, when each 

contingency analysed using the D.F. model at the same operating 
point. If weighing factotns for P.I. calculation using D.F.method 
is considc^red as unknown, the performance index for kth contingency 

a n b e w t' i 1 1 e ti as 



“ 1 *ki“i 

i 



(3.8) 


i = 1 , . . , N for voltage P.I. 

and 1 = 1 , . . , Ng for reactive power P.I. 

A. , is the coefficient and defined froin Eqns. (3.4) and (3.5) for 

K X 

known value of n as 


A 


ki 


i 

n 


V. 


- 


n 


^ .,1 im 
AV . 

1 


for voltage P . I 


(3.9) 


and 




1 

n 


QG.- 


AQG^-^ 


for reactive power P.I 


(3.10) 


Thus the function to he minimized is of the form of 
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F 


N 

c 



k = l 



A, , . w . 

ki I 



(3.11) 


Fof applying LSr, and minimizing error, it is essential that k > i 
in the above equation, which is generally the case in a power 
system network. 


In vector notation, the error vector can be defined as 


A.U' - b 


(3.12) 


where A is the rectangular matrix containing coefficients A, , and 

<jl* 

iJ is the weight vector and b vector contains elements b 

k ' 

Using LSE approach the optimum values of weights can be 
computed as 

■1 


U 




T - 
A b 


(3.13) 


Using aV)Ove equations, the optimum weights for both voltage 
performance indices Eqn. (3.4)] and reactive power 
performance Indices Eqn. (3.5)] can be computed. 


3. S SYSTEM STUDIES AND RESULTS 


Simulation studies were carried out on the 14 bus and 19 bus 
systems as described in Appencides A and B to test the potential 
of the proposed performance indices defined in Eqns . (3.4) and 
(3.5) for contingency selection (contingency ranking). All the 
contingency cases as considered for the two systems in the 
Chapter-2 were simulated. Results of the voltage contingency 
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analysis obtained by the exact method and distribution factors 
were utilized to compute the voltage performance indices assuming 
ail t.he wiMghtages to be unity. For computing the reactive power 
pet £ ormatice indices, the load flow was run without checking 
Q- limit oi generators and distribution factors were computed using 
these results as explained in Chapter~2. Contingency analysis 
results Wet e obtained using these sets of distribution factors and 

also by rutinitig exact laQ,d flow without considering generator 

Q *■ J 1. m 1 1 a , 

Dif ferent forma of voltage and reactive pOtSer performance 
indices wet e used considering the values of n = 2, 10 and 20 in 
Eqt»s t 1 - ^ ) and (3.5), 

The above three seta of voltage performance indices for 
IEEE -14 bus system using distribution factor’s results are 
presented in table 3.1(a) and those obtained from the exact load 
flow results are given in table 3.t(b). The, reactive power 

performance indices using distribution factors and exact load 
flow results are presented in tables 3.2(a) and . 3.2(b) 

respectively. Since exact load flow considering generator 
did not c'ouverge for line-1 outage, the voltage performance index 
(VPI) is not included in table 3.1. The relative severity index 
of various con t irrgenc i es are also Included in the above tables. 
It can be obser/'ved from these tables that with second order 
voltage performance index some of the severe contingency cases 
have been ranked lower than those in which there is no violation 
For IEEE~14 bus system the operating limit of 


of voltage. 
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wa^ considered as 0.95 to 1.05 p.u. whereas in the Indian 
.system, it: was taken as 0.9 p.u. to 1.1 p.u. This masking effect 
giadually t educes with the use of higher order performance indices 
and gets eliminated when exponent n = 20. However, in case of 
reactive power performance index (QPI) there is no significant 
change in ranking by using higher than second order performance 
ttidii'es, .hifliilar obs et’vat i one can also be seen in case of 19 bus 
Indian .system as is evident from tables 3.5Ca) i 3.5(b)and 3.6(a) 
and 3 . 6 ( b ) . 

In the above study all the welghtages (w , and w., , ) were 

vl QL 

f'onsideietl to be unity in both distribution factors method and 
exact load flow method. In fact one could have used user specified 
weightages Instead of equal weightages considered in the study. 
The values of performance indices computed with the distribution 
fai.*ti»ra and the exact load flow even with higher order exponent 
ate quite different and provide different rankings to the 
contingencies. In order to eliminate this misranking problem, the 
optimal values of weights to be used along with the distribution 
factors method were computed using least square error (LSF) method 
fot both VPl and QPI in case of IEEE-14 bus system only. The 
optimal weights SO obtained are given in tables 3.3. The 
performance index values obtained from exact load flow and used in 
LSE method were their earlier values obtained with all the weights 
cousldet ed as unity. The new values of voltage and reactive power 
perfoimance indices for the 14 bus system using distribution 
factors were computed and are compared with those obtained with 
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thf,> ‘ixai-t Load flow method as given in table 3.4. The new 
performance indices are quite close to those obtained using exact 
load flow, thus elimitxating the misranlcing effect. 

3.6 CONCLUSIONS 

In this chapter the higher order performance indices have 
been proposed and tested on two sample systems. Optimal weights 
Were obtained using LSE method and were used’ for contingency 
ranking Ttie* system study reveals that 

(. J .) The use of higher order exponent in voltage performance index 
el I mi «i.at e« the masking effect. However, it does not Improve 
the lesults of reactive power performance index. In voltage 
pel ftn mance index an exponent of 20 or more is recommended to 
ovetfome the masking effect. 

(2) tilth the optimal weights obtained from LSE method, the 
performance indices computed by the distribution factors 
matches with those obtained by the exact load flow method 
thus mlnlmiaing the mlsranking effect. 
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TABLE-3. la 


i.UMH.ARISOM OF VP la OBTAINED BY FACTORS METHOD 
UlTH VARIOUS EXPONENTS 


Out 

L 1 n# n 

VPI C‘ 

VPI 10 

VPI 20 

Bus Voltage 
St at us 

Lint No 

1 C 


.226510 

. 1 04743 

— 

2 S 

? 1 -) 0 00 1 0 ,1 ) 

2.778026(03) 

14.833389(03) 

V 

3 4 

t 6 3811(04) 

1 ,7858-41 ( 04) 

6.835161(04) 

V 

4 2 

\ B 18 9(06) 

.212892(06) 

.099986(07) 

V 

B i 

i'. 'i 1 0^'. ) 

687767(05) 

1 . 021348(05) 

V 

4 1 

U / 9 / ( 1 ) 

029307(13) 

.001947(14) 

NV 

7 1 

L / ( 0 b y 1 t .1 ) 

028740(14) 

. 002040(13) 

NV 

S v! 

0,.. ^80 1 ( Ob ) 

1 12868(08) 

. 045080(08) 

NV 

9 

u 0 4 ^ 0 c 0 9 ) 

. 068496(1 1 ) 

.010516(11) 

NV 

1 0 i 

‘ . *1 V 9 0 1 i t f ) 

022591(17) 

. 001142(17) 

NV 

i I 1 

bi)S‘.896 Ml) 

050743(12) 

. 007924( 12) 

NV 

1 2 1 

Ifa) 

, 021238(18) 

. 00096408) 

NV 

13 1 

HI i y j 0 1 1 0 i 

224596(07) 

. 232471 ( 06) 

V 

1 4 

‘ly Ob /b ( 2 1 ) 

009908(21 ) 

. 000254(21 ) 

NV 

15 1 

1 ?l u60(£0) 

012174(20) 

. 000349(20) 

NV 

1 6 1 

4 3 4S4 9( 1 9 ) 

. 019421(19) 

. 0006981 19) 

NV 

1 7 2 

0 34 1 sac 0?) 

.0776131 10) 

. 0132561 1 0) 

NV 

tS 1 

s 2 a S 6 S C 1 ? ) 

. 024652(15) 

. 001375(15) 

NV 

19 1 

£0206S( 18) 

, 023330(16) 

. 001262(16) 

NV 

20 t 

6 7c:9c‘0C 1 4 ) 

. 088927(09) 

.031188(09) 

NV 

Cen No . 

£ ? & 

Oufcrse ( 02 ) 

.101415^0(02) 

. 15848E+03{02) 

V 

3 12 9 1 227 (-( 0 1 ) 

. 150245+0(01 ) 

. 36853E+0S( 01 ) 

V 


No V Violating 

N V Non- V 1 0 1 at 1 ng 

NmnbiM givon in prantha«is refers to ranking. 

VPI c Voltage Performance Index with exponent 2 
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TABLE 3.1b 


tOhr mK 1 ivJN Of- VHIs OBTAINED BY RUNNING FULL AC LOAD FLOW WITH 

VARIOUS EXPONENTS 


Ou t .1 

i l ! > i 

4 ^' 

^ G i‘ f 1 


vr 1 2 

VPI to 

VPI 20 

Bus 

Voltage 
St at us 

1 J t If 

1 

H V 

i 

4 

vJ How h#* 

G V 0 3 4 { 3 ) 

not converged 
.27H0E+0I (S 

. 1 126SE + 02C6 ) 

V 

1 


I* 

/n9G3(3) 

.a930aE+03( 1 ) 

. 39303E + 07( 1 ) 

V 



1 

aoe548( 16) 

.59S94E-01 (17) 

.B9248E-0£( 17) 

NV 

e; 


c* 

,337c'£?( 10) 

.210342+00(11) 

. 1 066 1E + 00( 1 1 ) 

V 

It 



^l!SS00a(9) 

. 14606E+00( 12) 

. BSOiaE-OI ( 12) 

NV 

1 


s 

OaOSTTI 13) 

. 1 01 01E+00( 13) 

.22385E-01 ( 14) 

NV 

a 


4 

156054(4) 

. 19433E+02(3) 

.121622+04(3) 

V 

9 , 


1 

6£?647(20) 

.34638E-01 (20) 

.246l6E-0£(20) 

NV 

1 ii 



1 !■ 6308( t a) 

. 80009E-01 (14) 

. 15793E-01 (15) 

NV 

! 1 



S48T?S(6) 

.207115+01 (6) 

, 1254SE + 02(S) 

V 

1 i/ 


1 

I8( 1 9 ) 

76727E-01 (IS) 

.£3123E~01 (13) 

NV 

1 1 


It* 

114S42( 1 1 ) 

. 102S7E+02(4) 

. 52537E-f03(4) 

V 

1 4 


3 

1431 95(7) 

.40340E+00( 10) 

. 34086E + 00( 1 0 ) 

V 

1 


1 

.»96ia3( IS) 

.61 l30£-0t(16) 

. 1 0723E-01 ( 16 ) 

NV 

1 6 


1 

7t7426( 18) 

. 41 983E-01 (19) 

.S8415E-02( 18) 

NV 

1 7 



747381 (8) 

.47540E+00(9) 

.54846E+00(9) 

V 

1 a 


1 

7a3oaa( i 7) 

.43081E-01 (18) 

. 45938E~02( 19) 

NV 

1 ^ 


1 

sissai (ai ) 

.25451E-01 (21 ) 

. J5738£“02(£1 ) 

NV 

iiO' 


t 

0?03a9( 14) 

.807832+00(8) 

. 31501E + 01 (7) 

V 

1 1 

Ho 

I 

tt0(>4?c'(5) 

. 1 ia70E+01 (7) 

. 1 9539E + 0I (8) 

V 



8 

, 371559(1 ) 

.30097E+02(2) 

. 13856E + 04(2) 

V 


N o i f V Violating i 

Nv Non-Violating 

Number given in parenthesis refers to the ranking. 
VPI c.' Voltage Performance Index with exponent 2 



TABLE-3. ea 


COHPAR ISON 

OP OPIs OBTAINEP BY 

FACTORS METHOD 

UltH 

VARIOUS EXPONENTS 

Outage 

QPI a 

QPI to 


QPI £0 

L i n e / G e n 

«... .■» ... .... ..... 





tine No 

1 

S8t8SE+01 (3) 

.S9578E+03(3) 


.901 17E+06(3) 

£ 

33621E+01 (20) 

.261 1 6E+03(5) 


. 34058E+06 (5 ) 

3 

6 1 466E+01 (2) 

. 15464E+05(£) 


. 1 1 957E+1 0(2) 

4 

3491 lE+01 (17) 

.25833E+03(6) 


.33305E+06(6) 

S 

34H0E + O1 (19) 

.a66t7E+03(4) 


.3S384E+06(4) 

6 

3sa7-qe+oi (7) 

.85538E+03 


. 32S53E+06 

j 

aSSTSEt-OI (6) 

.a5538E+03 


. 32553E + 06 

a 

37796E+01 (5) 

.a5539E+03 


.325535+06 

9 

69812Et01 ( 1 ) 

.23507E+05( 1 ) 


.276295+1 0( 1 ) 

1 0 

3^95^Et-01 (16) 

.25539E+03 


.325535+06 

1 ! 

3S088E+01 (14) 

,25537E+03 


.325535+06 

1 a 

asaaaEt-oi ca ) 

.25539E+03 


. 32553E + 06 

1 3 

3S200E+01 (11) 

.25538E+03 


. 325535+06 

1 4 

408-16Et-01 (4) 

.25572E+03(7) 


. 325535+06 

IS 

351 31E+0t (13) 

.25538E+03 


. 325535+06 

1 6 

ISOBYEtOI ( 15) 

.25S37E+03 


. 325535+06 

1 T 

34790E-*-01 (18) 

.25661E+03 


. 32876E + 06 

1 8 

asai^E^'Ot (10) 

.25539E*f03 


. 325535+06 

1 9 

35a3lE-f01 (9) 

.255395+03 


.325535+06 

ao 

. 351 85E+01 ( la) 

.25538E+03 


.325535+06 

Gtrn No . 

a 

.a4997£+oi (aa) 

.a4871E+03 


.309275+06 

3 

. 31 546E*f-01 (81 ) 

.25454E+03 

r-' 

.323615+06 


Not®. QPI 2» Reactiv® Power Performance Index with exponent £. 
Number given in parenthesis refers to the ranking. 



6 


TABLE 3,Si& 

■'.k i Sou OF QFlb OBTAINED AFTER RUNNING FULL AC LOAU FLOW 
WITH VARIOUS EXPONENTS 


L. i I s e b •' 1 1 

0 li t o g e 

(IHl C 

QPI 10 

QPI SO 

Line f‘lv 

1 

b 0 V H 3 L * 0 1 (OS) 

. 3481 0E+04( 02) 

. 6(i498E+08( 02) 

c’ 

50S07F.fOl (03) 

. 5S241E+03( 04 ) 

. i4538E+0?(04) 

:s 

fO&oSEtOl (01 ) 

.20961E+05(0) ) 

.2f967E+1 0(01) 

4 

J 3 0 0 1 £ <■ 0 1 ( 1 0 ) 

. 521 S2E+03( 05 ) 

. 1 3S96E+0?( 05 ) 


4 1 4 0 } £ *' 0 1 ( 0 e ) 

. 91393E+03(03) 

.4i756E+07(03) 

6 

3‘".335Et 01 ( M ) 

. 321 91 E+03( 1 0) 

. Sl 725E + 06 (10) 

7 

3‘.:.;3S3E + 01 (17) 

.S6414E+03( 14) 

. 3^3E2E+06 (14) 

3 

3C I 7::SEtO 1 (13) 

.3E755E+03(08) 

.53550E+06( 03) 

9 

S V- 4 U 4 E 0 1 ( S S ) 

, 67852E+0t (£2) 

. ti916E + 03(22) 

1 0 

3E9SSEtO 1 ( 08 ) 

. 32646E+03( 09 ) 

.5il99E+06(09) 

1 1 

3 9 4 4 9 £ •♦ 0 U 0 7 ) 

.28184E+03( 13) 

•.39642E+06(13) 

1 u-' 

3 u il 8 £ * 0 1 ( 1 S ) 

.208042+03(20) 

.2i598E+06(20) 

\ 3 

3££84E 1 (It) 

.205152+03(21 ) 

.21 002E+06 (21 ) 

1 4 

4 HSI'lEi Ot ( 04) 

.23672E+03{ 18) 

.2f596E+06( 18) 

# j 5 

M9 3SE * 01 (19) 

.25255E+03( 17) 

.3I831E+06(17) 

1 6 

34 1 0 3£ * 0 1 ( SO ) 

22573E+03( 1 9 ) 

2S430E+06 (19) 

1 T 

4 1 SnGE i 0 1 ( OS ) 

. 49625E+03( 06 ) 

. 12242E + 07 ( 06 ) 

1 B 

3£497l£i'01 (16) 

. 28612E+03( 12) 

.46862E+06(12) 

1 9 

3SS39E1 01(18) 

. 25722E + 03( 1 6 ) 

.33023E+06(16) 

0 

3‘.S<;3c* 01(15) 

.29494E + 03( 1 1 ) 

. 43421E + 06( 1 1 ) 

e n N 0 

b 1 „3 „S ’/ £ ■»■ 0 1 (SI ) 

.26£51E+03( 15) 

.34377E+.06( 15) 

3 

1 8 1 0 u £ t 0 1 (09) 

. 3S596£ + 03< 07) 

.63I40E+06(07) 


Note QP I & - Reactive Performance Index obtained with 

exponents. 
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TABLE-3. 3a 


OPTIMAL WEIGHTS 


ASSOCIATED WITH 


BUS VOLTAGES 


Bus No 


? 

8 

9 

1 0 
t I 
1 £ 
I 3 
1 4 


Optimal Weight 


.24341768E+1 1 
-.806S9984E+08 
. 69220720E+06 
-.862215t3E+13 
-.87613909E+09 
.57939656E+08 
15597948E+05 
13392826E+05 
-.487771 07E+07 
.41619617E+08 
-. 127S8524E+08 
.47526336E+07 
.22186640E+04 


Note Optimal weight associated with bus-l(slack) is 

i no i gnif i cant as the voltage deviation is zero at 


OPTIHAL 
Q*™Source No 
I 

£ 

3 

4 

5 


TABLE~3.3b 


. 43942874E+00 
. 95915733EHHOO 
. 82600796E+00 
.49014612E+00 
.56535445E+0^ 


UEIGHTS ASSOCIATED WITH Q~S0URCES 

Optimal Weight 


slack . 



TABLE 3.4a 

COMPARISON Of VPIs OBTAINED BY FACTORS METHOD USING THE 
OPTIMAL WEIGHTS WITH THAT OF FULL AC LOAD FLOW 


L 1 / G e r'l 

out age 

Factors method 

AC load flow 

L i e No. 

1 

769 . 8601558927083 

not converged 


1 1 . 265-4247291502S 

1 1 .265376 

1 

930296 . 499956145 

393ii296.5 

4 

£ . 39 61 7680448 134£E“0£ 

5.9t4794000000000E-03 

S 

. 1 042288454507342 

.10661497 

B 

- . 123387170798395 

3.5011 94299999999E-02 

7 

-e 967237402279959E-02 

2.2384771 OOOOOOOOE-02 

8 

1216.206839859732 

1211.1917 

9 

-7 9791 1719787953 

£.4il574600000000E-03 

1 0 

S . 1 1 9544230896732E-02 

1 .579288200000000E-02 

1 1 

12 . 66541614815903 

12.^4927 

1 a 

-1.1 5279716961427 

2.3i23£4400000000E-0£ 

1 1 

SI 6 , 9864674322652 

525.36812 

i 4 

2381854488629256 

. 3408643 

1 S 

2 07592603899705 

1 . ot£326300000000E-0£ 

1 4 

772S866222321227 

5.841471 600000000E-03 

1 7 

1 00716034691231 

. 54846344 

1 8 

1 1 5277060416016 

4 .593761 700 OOOOOOE-03 

1 9 

2 . 31509199488985 

1 . 573764000000000E-03 

ao 

67 . 990316431517 

3.1501054 

C 1 11 N 0 

1 , 9S387060800567 

1 .95387063 

3 

1385.563001535832 

1385.563 


i 



TABLE 3.4b 

COMPARISON OF QPIs OBTAINED BY FACTORS METHOD USING THE 
OPTIMAL WEIGHTS WITH THAT OF FULL AC LOAD FLOW 


Lin©/G^n AC Lodd Flow Factors niethod 

outage 


Line No 

1 

S 81846920885382 

5, 6^3 

£ 

3 . 36214739474536 

5. 0fe07 

3 

6 . 1 4659428889309 

7. 0868 

4 

3.491 12064753459 

3.8601 

S 

3.41403665125746 

4 . 1^07 

6 

3.58737739914148 

3.5885 

7 

3.58780988724758 

3.5388 

8 

3.77959992436361 

3.6125 

■9 

6 . 98123891 1712 

£ . 6404 

1 0 

3.49541131742889 

3.8625 

1 1 

3.5087966387337 

3.9449 

1 8 

3.58319658549863 

3 . 63l 8 

t 3 

3 , 5200317172085 

3.6S84 

1 4 

4 . 08460713161885 

4.3il3 

1 S 

3.51308490859184 

3.4635 

t 6 

3 50878584963889 

3 .4508 

1 7 

3 .47901835380895 

4 . 1 S 88 

1 8 

3 .S818170988574 

3.5497 

1 9 

3 .58305483743849 

3.5239 

ao 

3.51846493534591 

3.5563 

Gen No 

2 

2.49966072210234 

3.1337 

3 

3. 15461951397959 

3.81 



TABLE-3 .5a 
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COMPARISON 

OF VPI* OBTAINED BY FACTORS METHOD 
WITH VARIOUS EXPONENTS 


vu'i 

VPI 2 

VPI 10 

VPI 20 

Bus Voltag# 

t. ina/Can 




Status 

1 11 # No . 





1 

. 32S87E+0E 

. 13521E+05 

. 1T651E + 0f 


e 

. 1 f 1 62E4-(|2 

. 1 0612E + 04 

, I323EE+07 


3 

. 1 5240E+02 

. 34023E+03 

. 143f4E-*.06 


4 

. S6576E+01 

. 74292E+01 

. 1 0309E+03 


s 

.22274Et01 ( 14) 

.5181£E-01 ( 14) 

.584722-02(14) 

NV 

6 

. 1 9992E+03( 03) 

. 88402E+08( 03) 

,80166E+16(02) 

V 

T 

76407E+03( 01 ) 

.73101E+1 1(01) 

.58E54E+22(01 ) 

V 

8 

20282E+01 (IS) 

. 31247E-01 (15) 

, 20670E-02( 15) 

NV 

9 

. 244TIE'f0! C 13) 

.81771E-01 ( 13) 

. 15159E-01 ( 13) 

• NV 

1 0 

. 5E046E+02( OS) 

. 1297SE + 06( OS) 

. 16585E+1 1 ( 05) 

V 

1 1 

.13211 E + 02( 1 0 ) 

. 19519E + 03( 10) 

. 6741 1E+05( 1 0) 

NV 

1 e 

. T3T41 04) 

. 67069E + 06( 04) 

. 4271 SE-^-l 2C 04 ) 

V 

13 

3S372e+02{ 07) 

. 1T560E1-05C OT) 

. 43T65E4'09( 07) 

V 

1 4 

. 20577E+03( 02) 

. 938552^08(02} 

. 76386E+1 6( 03) 

V 

15 

. 1 06242+03 

. 4l502Ei*08 

. 24570E+1 6 


1 5 

. 20497E+02{ 08) 

. 94185E+04{ 08) 

.21 044E+09( 08) 

V 

IT 

, 35623E+02( 06) 

.E6019E+OS( 06) 

, 16015E4*10(06) 

V 

18 

. 1 9439E+02( 09) 

. 16676E + 04(09) 

. 66885E+07( 09) 

V 

If 

.24275E+03 

. 13955E+09 

.28274E+17 


£0 

. 1 6851E+03 

.36289E+08 

. 13205E + 1 6 


£1 

. 20307E+0) 

.29660E-01 

. 17461E-02 


££ 

1 1 163E+03 

. S0967E+07 

.33818E4-14 


23 

685002+02 

. 681 13E+06 

.63651E-H2 


24 

. 35520E+01 

. 67213E+00 

. 1 061 OE+01 


No . 





£ 

.287I9E+02 

.23542E+04 

.57747E+07 


3 

. 9E881E+01 (11) 

.43E35E+02( 1 1 

.32297E+04( 1 1 ) 

NV 

4 

. 43062E+01 (12) 

. 13896E+01 (12) 

.43890E+01 ( 12) 

NV 

c # V 

Violating 





NV •• Non-Violating 

Numbar givan in pranthaais rafars to ranking. 

VPI Voltaga Parformanca Index with exponent £ 
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TABLE 3.5b 

COMPARISON OF VPIs OBTAINED BY RUNNING FULL AC LOAD FLOW WITH 

VARIOUS EXPONENTS 


VPI S VPI 10 VPI SO Bt 

Vc 

St 


.in* No. 


5 

.5707£e4E+01 (09) 

.a5S040aE+0i (i o) 

. 1 0H7489E+02C 1 0) 

h 

§ 

. 1 943457E+04t OS) 

.S85£177E+13(0S) 

.6271 t53E+a5(02) 


7 

. 1E91 798E+03( 04) 

.E6£7690E+07(04) 

.43T1T30E+13C04> 


8 

. 1388364E+01 (15) 

.597S43EE-0S( 14) 

.8T96558E-04(13> 

N 

9 

.S4S9371E+01 (IS) 

.S731S5EE-01 (IS) 

.801S052E-02( 12) 

N 

1 0 

.S478646E+0a(06) 

. 13130S7E+05(06) 

.4138272E+09C06) 


1 1 

. 1 9830S8E+01 (13) 

.9401519E-0S(13) 

,7152308E-04C1 4) 

N 

IE 

559ES39E+0S( 05) 

. 1S530taE+06(05) 

. 1832864E+11 COS) 


13 

.S360941E+04(01 ) 

.6149688E+15{01 ) 

.3315693E+30( 01 ) 


14 

161 0369E+04( 03) 

. 1 13£791E+13(03) 

. 95T3743E+24C 03) 


1 6 

.56S7973E+01 (10) 

. 1 18S84SE+0S( 09) 

.359S399E+a3C09) 

N 

IT 

.866E158E+01 (08) 

.18880995+03(08) 

.2060990E'»'05( 08) 

N 

ta 

34S6S67E+01 (11) 

.3704335E+00( 1 1 ) 

.5060231E+O0C 1 1 ) 

N 

No . 
3 

14eES96E+01 ( 14) 

.367E53£E-0S( 15) 

. 1 189609E-04( 15) 


4 

, SI 80394E+0E( 07) 

. SS66990E+04( 07) 

. 7452152E+07( 07) 




TAeLE-3.6« 


COMPARISON OF QPI, OBTAINED BY 


factors method 


yui 

QPI 2 

1. irf#/C»n 




i ri» No . 

1 

48088E+02 

£ 

. 22395E+01 

3 

20405E+02 

4 

536265+01 

5 

14874E+01 ( 09) 

B 

. 491 30E+03 

7 

1 99305 + 04 

a 

1 62855+01 ( 08) 

9 

• 1 46405 + 01 ( 10) 

0 

971995+02(02) 

1 

18324E+02( 06) 

£ 

144405+03 

3 

715855+02 

4 

50182E+03 

S 

552475+03 

. 6 

83787E+02(03) 

1 7 

1 1 256E+03 ( 0 1 ) 

i 8 

517705+02(05) 

f 9 

77966E+03 


37526E+03 

'I 

145455+01 


313325+03 


.216725+03 

£4 

328765+01 

C*I| No . 

s 

. 794905+02(04) 

3 

. 1 70705+02 

4 

17076E+02(07) 


QPI to 


•70779E+09 
. 1Ee53E+01 
.78818E+07 
• 30364E+04 
.a5477E+00 
. 90251E+14 
■ 1 0OS6E+1 6 
■ H5532E+00 ( 08 ) 
.25477E+00 
.25750E+I 1 (02) 
.437S9E+07(06) 
• 1 91 02E+12 
. S4297E+1 0 
. 1 0037E+15 
. T6E5SE+15 
.121 08E+1 1(03) 

. 541 79E+1 1(01) 

. 1 0322E+1 0(05) 

. 91336E+15 
. 23355E+14 
. e5477E+00 
. 94402E+13 
. 14794E+13 
. 66995E+02 

. 95533E+1 0( 04) 
.34638E+07 
. 34638E+07( 07) 


WITH VARIOUS EXPONENTS 
QPI 20 


. 25048E+1 9 
.48179E+01 
.31061E+15 
. 46091 E + 08 
• 1 0821 E+00 
.40726E+29 
. 50263E+35 
. 1 0821 E+00 
. 1 0821 E + OO 
.331522+22(02) 
.95742E+14(06) 
. 18245E+24 
. I4741E+21 
. 50372E+29 
. 1321 OE+30 
. 73306E+21 ( 03) 

. 1 4677E+23( 0 1 ) 
.532755+19(05) 
.4171 lE+31 
. 27274E+28 
. 1 0821E+00 
. 44559E+27 
. 1 0944E+26 
.22271E+05 

. 45632E+21 ( 04) 
.59989E+14 
. 59989E + 1 4(07) 


TABLE 3.6b 

COMPARISON OF QPI» OBTAINED AFTER RUNNING FULL AC LOAD FLOW 

WITH VARIOUS EXPONENTS 


Outag* 

.in* No. 
5 
8 
9 

1 0 
1 1 
1 6 
17 
16 

. wn No . 
E 
4 


QPI 2 


QPI 10 


QPI 20 


151 04E+0E(03) 
21 28eE+01 ( 09) 
19716E+01 (10) 
2T178E+02(02) 
491S2E+01 (07) 
56092E+01 ( 06) 
6777SE+01 ( 04) 
43534E+01 ( 08) 


14e38E+03(01 ) 
6015EC+01 (05) 


. 156455+07(03) 
. 726815 + 00( 09) 
. 38620E+00( 10) 
. 368685+08(02) 
. 16755E + 04( 07) 
.420835+04(06) 
. 144465+05(04) 
. 701375+03(08) 


. 18140E+1E( 01 ) 
. 1 05245+05(05) 


. 122395 + 14(03) 
. 124385 + 01 ( 09) 
.20632E+00( 10) 
. 679645+1 6(02) 
. 14033E+08( 07) 
.88538E+08( 06) 
. 1 04335+1 0( 04) 
. 245805+07(08) 


. 1 64525+24(01 ) 
.553745+09(05) 


CHAPTER - 4- 


CONCLUSIONS 


Power system security monitorifig and analysis forms ai 
integral part of the modern energy management system (EMS) but iti 
on”"line implementation poses a major challenge to the power aystei 
engiiieers , Recently the voltage security problem has found mor« 
attention by researchers. Contingency selection is performed t< 
liiiniml^se the, total number o£ cont ingencies to be analyzed by full 
At load flow and hence reducing the total time of security 
analysis. From the available literature on voltage contingency 
selection, it is observed that further improvements are needed ir 
the available methods, 

a) to adopt better performance indices which can truly reflect 
the relative severity of the contingency cases; thus 
minimizing the misranking and masking effects, inherent with 
existing second order performance indices 

b) to explore faster and more accurate models for predicting the 
post outage voltages and reactive powers for all cont ingenci es 

ht£ove ranking. j 

Keeping the above two objectives in mind, in the present 
thesis a new set of voltage and reactive power distribution 
factors have been suggested, which can be directly obtained from a| 
baae ease load flow results using sensitivity properties of N.R.| 



load flow Jacobiaiis. Higher order voltage and reactive power 
p e r f o t mane e indices have been tried out and a new approach for 
finditig optimal weights used in the performance indices, based on 
least square error method, has been suggested. System studies to 
test the effectiveness of the new distribution factors and the 
performance indices, were conducted on IEEE-14 bus and a 19-bus 
Indian system (representing UP State 400 kV transmission network). 
Following main conclusions can be drawn from the works reported, in 
this 

1) The proposed voltage distribution factors predict the bus 
voltages with sufficient accuracy. However, the proposed 
reactive power distribution factors are not equally accurate 
in predicting the reactive power outputs of sources. 

fKe.i'!!’ 

2) The calculation of distribution factors and - it - s updating is 
extremely fast as it utilizes only the available base load 
flow results. For alight variation in loading (upto ± 5%) the 
factors need not be updated and the one already computed can 
be used for predicting the bus voltages. This further reduces 
the computational time for security analysis. 

3) Use of higher order voltage performance indices eliminates 
the problem of masking effet. An exponent of 20 or more| 
should be used in the voltage performance index. However, use; 
of higher, exponents In reactive power performance index do not 
significantly affect the relative ranking of contingencies. 
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4) Use of the optimal weights obtained through the proposed LSE 
method in the voltage as well as reactive power Performance 
indices brings their values computed by the distribution 
factors method and the exact methods quite close. Thus, 
minimizes the misranking effect. 

5) For voltage contingency selection, it is recommended to , 

^ US ^ 

the proposed voltage performance index algorithm alongwith the 
suggested voltage distribution factors for post outage voltage 
prediction. This model is found to be superior to the 
reactive power performance index baaed model utilizing the 
reactive power distribution factors. The voltage contingency 
selection algorithm, so obtained, will be extremely fast and 

can be effectively used for real time voltage security 
analysis. 

As a consequence of the works carried out in this thesis 
further research efforts can be made in the following direction 

il The distribution factors suggested in this thesis hSiVe been 

derived for only a single outage of line/transformer 

generator. The distribution factors, on similar lines r. 

» can 

be explored for multiple contingencies. 

ill The reactive power baaed ranking has not been found effective 
in the present work. Further investigations are required to 
evolve better models. A hybrid performance index considering 
both bus voltages and reactive powers can also be tried out 
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iii) The studies have been conducted on only 14-bus and 19-bus 
systems. It can be tried out on larger practical systems 
containing several hundred buses and lines. 

iv) Recently neural networks have been popularly used for various 

power system studies. Being non-algorithmic in nature and 
extremely fast, it is ideally suited for security analysis 
problems. The neural network baaed models for the 

contingency selection and analysis can be tried out. 
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APPENDIX -A 

DATA FOR IEEE- 14 BUS SYSTEM ( 100 MVA BASE ) 
The IEEE- 14 bus network is shown in fig-A. 

Table-A. 1 

(Generator Bus Data) 


Bus no Scheduled real Specified voltage Reactive generation limit 
power generation magnitude (p.u.) 



Pg(MW) 

Vspe 

Max(MVAE) 

Min(MVAE) 

1 


1.060 

-- 


2 

40.0 

1.045 

0.5 

-0.4 

3 

20.0 

1.070 

0.24 

-0.06 

4 

0.0 

1.010 

0.4 

-0.0 


5 0.0 1.090 0.24 -0.06 


Real 

and reactive 

power load at generator buses are 

Table-A. 2 

(Load Bus Data) 

zero. 

Bus no 

Load 

External 

shunt 

susceptance 





(p.u. 

) 


Real (MW) 

Reactive ( MV AR) 



6 

0.0 

0.0 


0.0 


7 

29.5 

16.6 


0.19 


8 

7.6 

1.6 


0.0 


9 

47.8 

3.9 


0.0 


10 

9.0 

5.8 


0.0 


11 

3.5 

1.8 


0.0 


12 

6.1 

1.6 


0.0 


13 

13.5 

5.8 


0.0 


14 

14.9 

5.0 


0.0 
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Table-A. 3 
{Transformer Data) 


Bus no 

From 

To 

Series 

impedance 

Tap setting 


bus 

bus 







R 

X 





(P-u. ) 

(p.u. ) 


14 

8 ■ 

3 

0.0 

. 25202 

0.962 

15 

9 

6 

0.0 

.20912 

0.978 

16 

9 

7 

0.0 

.55618 

0.969 




Table- 

A. 4 





(Line Data) 


Bus no 

From 

To 

Series 

impedance 

Shunt susceptance 


bus 

bus 



B(p.u.) 




R 

X 





(p.u. ) 

(p.u. ) 


1 

1 

2 

.01938 

.05917 

.0264 

2 

1 

8 

.05403 

.22304 

.0246 

3 

2 

4 

.04699 

. 19797 

.0219 

4 

2 

8 

.05695 

. 17388 

.017 

5 

2 

9 

.05811 

. 17632 

.0187 

6 

3 

11 

.09498 

. 19890 

.0000 

7 

3 

12 

.12291 

. 25581 

.0000 

8 

3 

13 

.06615 

. 13027 

.0000 

9 

4 

9 

.06701 

.17103 

.0173 

10 

6 

5 

.00000 

.17615 

.0000 

11 

6 

7 

.00000 

. 11001 

.0000 

12 

7 

10 

.03181 

.08450 

.0000 

13 

7 

14 

. 12711 

.27038 

.0000 

17 

9 

8 

.01335 

.04211 

.00640 

18 

10 

11 

.08205 

. 19207 

.00000 

19 

2 

13 

.22092 

. 19988 

.00000 

20 

13 

14 

. 17093 

.34802 

.00000 




(a) 

H uA I E ££-14 BUSSYSTEM 

( Buses renumbered) 

O'ElINE, diagram 

3- WINDING URANSFOPMER EQUIVALENT 





GENERATOR 


SYNCHRONOUS 

CONDENSER 

SHUNT REXTOR/ 
CAPA Cr TOR 

transformer 




APPENDIX -B 


DATA FOR UPSEB 400KV-19 BUS SYSTEM ( 100 MVA BASE ) 
The UPSEB 400KV~19 bus network is shown in fig-B. 


Table-B. 1 

(Generator Bus Data) 


Bus no Scheduled real Specified voltage Reactive generation limit 
power generation magnitude (p.u.) 



Pg(MW) 

Vspe 

Max(MVAR) 

Min(M 

1 


1.03 


— 

2 

1600.0 

1.03 

144.0 

0.0 

3 

900.0 

1.03 

80.0 

0.0 

4 

550.0 

1.03 

50.0 

0.0 

Real 

and reactive 

power load at generator 

buses are 

zero . 


Table-B. 2 
(Load Bus Data) 


us 

uo 

Load 

External shunt reactance 


- — -■ 

..... .... ... _ 

(p.u. ) 


Real (MW) 

Reac t ive ( MVAR ) 


5 

0.0 

0.0 

.494395 

6 

191.0 

16.0 

1.10250 

7 

1000.0 

0.0 

1.75000 

8 

0.0 

0.0 

2.20500 

9 

135.0 

78.0 

2.20500 

10 

236.0 

45.0 

2.20500 

1 i 

360.0 

11.0 

.676380 

12 

337 , 0 

59.0 

2.20500 

1 3 

90.0 

20.0 

0.00000 

14 

520 . 0 

84.0 

2.20500 

15 

387.0 

64.0 

1.10250 

16 

288.0 

0.0 

0.40384 

17 

477.0 

00.0 

2.20500 

18 

"151.0 

69.0 

0.00000 

19 

-263.0 

00.0 

0.00000 



Table~B. 3 
(Transformer Data) 


Bus no 

From 

bus 

To 

bus 

Series 

impedance 

Tap setting 


R 

(p.u. ) 

X 

(P-u. ) 


1 

c. 

3 

4 

5 

6 

7 

8 

2 

1 

3 

4 

.00016 

.00073 

.00030 

.00049 

.00591 

.01460 

.01199 

.01943 

1.000 

1.000 

1.000 

1.000 


Table-B.4 
(Line Data) 


Bvis no Fi'om To 
bus bus 


5 

5 

7 

6 

7 

16 

7 

16 

5 

8 

8 

6 

9 

6 

8 

10 

B 

9 

li 

10 

9 

12 

6 

10 

13 

10 

■ 11 

14 

11 

5 

15 

11 

12 

16 

13 

12 

17 

13 

14 

18 

14 

12 

19 

15 

14 


Series impedance 


R 

X 

(p.u, ) 

(p.u. 

.00031 

.00310 

.00918 

.09306 

.00439 

.04464 

.00031 

.00310 

.00051 

.00517 

.00479 

.04880 

.00254 

.02584 

.00468 

.04770 

.00294 

.02997 

.00823 

.08375 

.00650 

.06617 

.00325 

.03307 

.00370 

.03762 

.00260 

.02646 

.00806 

.08169 


Shunt suaceptance 
B(p.u.) 


.04056 

1.21680 

2.34148 

.04056 

.06760 

.63614 

.33798 

.62450 

.39206 

1.09503 

.86521 

.43264 

.48870 

.34610 

1.06808 




